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Frontispiece 

The  composite  scene  above  pictorially  illustrates  many  of  the  radiometric  relationships  from  which  the 
implications  discussed  in  this  report  and  derived.  The  photographic  portion  illustrates  a  joint 
airbome/groundbased  experimental  measurements  program  conducted  in  the  vicinity  of  Mt.  Ranier, 
Washington  by  the  Visibility  Laboratory  of  the  University  of  California,  San  Diego,  in  cooperation  with 
and  under  the  sponsorship  of  the  Air  Force  Geophysics  Laboratory.  The  cartoon  overlay  illustrates  the 
influence  of  the  direct  and  diffuse  light  fields  upon  surface  targets  and  the  atmospheric  path  between  the 
surface  and  the  airborne  observer.  Within  the  visible  spectrum,  these  atmospheric  processes  are  primarily 
scattering  in  nature.  Whereas,  within  the  infrared  regions,  the  processes  must  be  expanded  to  include  the 
effects  of  absorption  and  emission.  A  discussion  of  these  combined  effects  and  their  interrelationships  is 
the  primary  content  of  this  report. 
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SUMMARY 


Thii  report  Is  a  part  of  the  Visibility  Laboratory's  Project  OPAQUE  eifort  and  was  prepared  under  AFGL  Contract 
F19628-82>C*0060.  It  describes  a  review  of  the  theoretical  constructs  for  a  rr^lation  medium  which  scatters,  absorbs  and 
emits  radiation  with  special  emphasis  on  applications  to  the  visible  and  the  infrared. 

The  equation  of  transfer  for  radiance,  irradiance  and  scalar  irradianr*  • .  fully  defined  and  then  integrated.  If  the  4ir 
radiance  distribution  is  measured  at  several  altitudes,  the  radiance  arrays  can  be  tested  for  consistency  Bnd,  if  consistent,  a 
measure  of  absorption  obtained. 

Some  of  the  special  characteristics  of  the  solar  almucantar  are  delineated.  Methods  are  presented  for  using  the  solar 
almucantar  sky  radiance  to  obtain  radiance  transmittance,  test  optical  stability,  and  test  for  spurious  sun  reflectance  in  the 
sky  radiance  photometer. 

Measurements  made  with  the  Visibility  Laboratory  integrating  nepheiometer  are  reviewed  and  compared  to  the 
Barteneva  (1960)  catalog  of  photopic  volume  scattering  functions.  An  average  volume  scattering  function  is  derived  and 
methods;  developed  for  obtaining  the  single  scattering  albedo  from  horizon  sky  radiance,  for  obtaining  the  scattering 
transmittance  from  sky  radiance  ratios,  and  for  obtaining  the  aerosol  optical  thickness  from  solar  almucantar  radiances. 
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IMPLICATIONS  OF  THE  EQUATION  OF  TRANSFER  WITHIN  THE 
VISIBLE  AND  INFRARED  SPECTRUM 

Jacqueline  I.  Gordon 


1.  INTRODUCTION 

In  the  increasingly  sophisticated  world  or  electro- 
optical  detection,  search,  and  guidance,  the  requirement 
for  establishing  and  predicting  atmospheric  influences  on 
system  performance  continues  to  be  a  primary  operational 
necessity.  It  is  In  support  of  this  requirement  that  the 
Visibility  Laboratory  in  cooperation  with,  and  under  the 
sponsorship  of  the  Air  Force  Geophysics  Laboratory  has 
maintained  an  extensive  program  of  airborne  optical  and 
meteorological  measurements.  In  recent  years  this  pro¬ 
gram  has  been  conducted  as  an  independent  but  coopera¬ 
tive  effort  (Johnson  ei  al .  (1979))  in  conjunction  with  the 
NATO  program  OPAQUE  (Optical  Atmospheric  Quanti¬ 
ties  in  Europe),  Fenn  (1978).  During  the  two  year  inter¬ 
val  spanning  the  years  1977  and  1978.  over  80  missions 
were  flown  documenting  the  vertical  structure  of  the  visi¬ 
ble  spectrum  total  volume  scattering  coefficient  in  the 
lower  troposphere  as  well  as  the  4rr  radiance  distribution 
of  the  ambient  light  field  at  several  altitudes.  Since  a 
thorough  awareness  of  the  vertical  structure  of  volume 
scattering  coefficient  Is  essential  to  the  prediction  of 
atmospheric  influences  on  contrast  transmittance  through 
this  regime,  these  data  have  been  presented  In  a  series  of 
technical  reports,  the  two  most  representative  of  which  are 
entitled  "Airborne  Measurements  of  Atmospheric  Volume 
Scattering  Coefficients  In  Northern  Europe,  Summer 
1978,"  Johnson  and  Gordon  (1980),  and  "An  Analysis  of 
Natural  Variations  in  European  Sky  and  Terrain  Radiance 
Measurements",  Johnson  and  Hering  (1981b). 

The  optimum  use  of  the  experimental  data 
presented  in  reports  such  as  those  referenced  above  is 
surely  to  establish  the  baseline  assessment  of  those  optical 
characteristics  most  influencing  slant  path  contrast 
transmittance,  and  to  develop  from  these  assessments 
realistic  predictive  models.  An  initial  effort  in  this  model 
development,  using  both  surface  and  profile  data  from  the 
OPAQUE  program  Is  discussed  in  Johnson  ei  al.  (1979), 
and  the  further  application  of  these  data  to  contrast 
transmittance  modelling  is  illustrated  by  Hering  (1981a). 

In  order  to  extend  the  methodologies  developed  for 
visible  spectrum  modelling  into  the  infrared,  it  seems 
appropriate  to  review  the  theoretical  constructs  for  a  radia¬ 
tion  medium  which  scatters,  absorbs  and  emits  radiation. 
These  constructs  are  being  used  to  evaluate  the  data  from 
the  atmospheric  optical  measurement  program  and 
develop  models.  In  addition  we  will  review  a  few  of  the 
data  sets  which  have  specific  modeling  Implications. 


The  following  section  will  deal  with  the  basic  mono¬ 
chromatic  equation  of  transfer  us  it  relates  to  radiance, 
irradiance  and  scalar  irradiance.  The  concepts  will  be  fully 
defined  and  related  where  possible  to  the  theoretical  work 
of  Duntley  etal.  (1957),  Preisendorfer  (1976)  and  Chan¬ 
drasekhar  (1960).  A  table  of  notationul  equivalencies  for 
the  Chandrasekhar  (1960)  notation  and  the  Visibility 
Laboratory  notation  is  given  in  Appendix  B. 

The  third  section  will  delineate  some  of  the  special 
features  of  the  solar  ulmucantar. 

The  fourth  section  will  review  the  measurements 
made  with  the  Visibility  Laboratory  integrating  nephelom- 
eter.  These  data  will  be  compared  to  the  Barteneva  (I960) 
catalog  of  photopic  volume  scattering  functions  and  an 
average  volume  scattering  function  derived.  Some  of  the 
implications  of  the  average  volume  scattering  function  will 
then  be  explored. 


2.  EQUATION  OF  TRANSFER 
2.1  Radiance 

The  basic  equations  developed  herein  are  mono¬ 
chromatic.  Applicability  to  broadband  sensors  will  be  dis¬ 
cussed  where  appropriate. 

The  most  basic  equation  in  radiation  theory  is  the 
equation  of  transfer  (Eq.  10  Duntley  ei  al.  (1957), 
Eq.  (3)  Sec.  3.15  Preisendorfer  (1976),  and  Eq.  (46) 
Chandrasekhar  (I960)) 


dL(z$.d>)l dr  -  -u(z)  Ll:,9,d>)  +  L, (*,#,</>).  (2.1) 

This  equation  relates  the  incremental  change  in  radiance 
dZ.(z,0,0)  at  altitude  z  in  direction  zenith  angle  0  and 
azimuth  0  over  the  incremental  path  length  dr  to  the 
attenuation  coefficient  <*(z),  the  radiance  L(z,6,d> )  and 
the  path  function  Z,,(z,0,0).  The  first  term  is  the  loss 
term  and  the  second  is  the  gain  term. 

The  attenuation  or  extinction  coefficient  a(z)  is 
equivalent  to  the  Chandrasekhar  (1960)  mass  attenuation 
coefficient  k  times  the  density  p.  The  attenuation 
coefficient  is  the  sum  of  the  total  scattering  coefficient 
s(z)  plus  the  absorption  coefficient  a(z ) 

or(z)  -  .<(*)  +  a(z).  (2.2) 
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The  path  function  £„(z,fl,<£)  is  equivalent  to  the 
Chandrasekhar  (I960)  mass  emission  coefficient  j  times 
the  density  p.  The  path  function  is  the  sum  of  a  scattered 
component  and  an  emitted  component  which  is  related  to 
the  absorption 

£.(z,», «)  -  L,s (*.«,<*)  +  £,fl(z).  (2.3) 

The  emitted  component  due  to  absorption  is  isotropic  and 
hence  is  shown  without  direction  modifiers. 

SsaMini 

The  scattering  coefficient  has  two  components.  Ray- 
leigh  or  molecular  scattering  is  highly  wavelength 
dependent  being  proportional  to  the  inverse  fourth  power 
of  the  wavelength  X~*.  Mie  scattering  us  is  less 
wavelength  dependent.  Both  tend  to  be  smooth  continu¬ 
ous  functions  with  wavelength. 

The  path  function  due  to  scattering  is  the  integral  of 
the  incoming  radiance  in  ail  4rr  directions 

(including  tne  sun  whem  appropriate)  times  the  volume 
scattering  function  <r(;,j3') 

<r(i.0')dn.  (2.4) 

The  scattering  angle  p'  is  equivalent  to  the  Chandrasekhar 
(1960)  angle  ©.  Ills  a  function  of  the  incoming  radiance 
direction  angles  & ,fi’  and  the  sensor  direction  angles  0,<£ 
as  follows 


coj/3'  -  sin®  sin®1  cos(ot'-<ft)  +  cosfl  cosfl' .  (2.5) 

The  integral  of  the  volume  scattering  function  over  An  Is 
the  total  scattering  coefficient 

s(z)  -  J<r(i,p)dn  .  (2.6) 


The  volume  scattering  function  is  equivalent  to  the  Chan¬ 
drasekhar  (1960)  phase  function  p{ cos0)  times  the 
attenuation  coefficient  a  divided  by  47r, 

Absorption  and  Emission 

Atoms  and  molecules  in  the  gas  phase  absorb  and 
emit  in  line  and  band  spectra.  The  atomic  spectra  are  line 
spectra  and  tend  to  be  at  the  shorter  wavelengths.  The 
band  spectra  are  molecular  and  tend  to  be  at  the  longer 
wavelengths.  The  contlnua  are  essentially  part  of  the  band 
spectra,  weak  but  broad  spectrally. 

Molecules  in  the  liquid  or  solid  phase  cannot  emit 
line  or  band  spectra  but  can  only  emit  in  a  continuous 
spectrum,  tls  distribution  of  which  is  determined  by  the 
ambient  temperature,  in  other  words  like  a  black  or  grey 


body  (incomplete  radiator).  This  implies  that  the  absorp 
tion  spectrum  for  liquids  and  solids  is  also  continuous 
spectrally. 

An  aerosol  is  defined  as  a  mixture  of  gas  and  small 
solid  or  liquid  particles.  Although  the  solid  particulates 
and  water  droplets  with  dissolved  particulates  absorb  and 
emit  like  a  grey  body  with  a  continuous  absorption  spec¬ 
trum,  the  gas  or  air  absorbs  and  emits  in  line  and  band 
spectra.  Thus  the  aerosol  absorbs  and  emits  like  a  spectral 
or  colored  body. 

Emission  Mechanisms.  The  principal  emission 
mechanisms  above  the  mesosphere  are  electrolumines¬ 
cence  and  chemiluminescence  resulting  in  line  and  band 
spectra.  These  emissions  are  called  airglow.  These  are 
important  in  the  visible  part  of  the  spectrum  at  twilight 
and  night  but  will  not  be  dealt  with  herein. 

The  principle  emission  mechanism  in  the  tropo¬ 
sphere  and  in  the  atmosphere  at  or  below  the  mesosphere 
is  temperature  radiation  which  is  photon  emission  caused 
by  atomic  or  molecular  collision.  The  atmosphere  is 
assumed  to  be  in  local  thermodynamic  equilibrium,  hence 
RirkhofTs  law  applies.  The  path  function  due  to  emission 
is  thus  [Eq.  (38)  Chandrasekhar  (I960)) 

£v(z)  -fl(rU(x,r)  (2.7) 

where  Z.(X,D  is  the  blackbody  radiance  at  wavelength  X 
and  T  is  the  temperature  in  degrees  Kelvin.  The  path 
function  due  to  emittance  at  and  below  the  mesosphere  Is 
highly  wavelength  dependent  and  should  be  dealt  with 
monochromatically. 

A  black  body  in  thermodynamic  equilibrium  absorbs 
and  emits  as  a  continuous  function  of  wavelength  and 
temperature  according  to  the  classical  equation  [Wolfe 
(1978)  Table  1.7] 

L(X.r)-r,/[ir\3(f*-l)]  (2.8) 

where 

C,  -  2 it  he 1  -  3.741382£-16  (2.9) 

and 

Jf-cj/(xn.  (2.10) 

The  form  3.74 E- 16  is  an  alternate  format  for  3.74xl0~16. 
This  computer  form  is  used  throughout  this  report.  The  c 
is  the  speed  of  light,  h  is  the  Planck  constant  and 

c,  -  ch/k  -  1.4387B6£-2  mK  (2. 1 1) 

where  k  is  the  Boltzmann  constant.  The  constants  are 
from  Driscoll  and  Vaughn  (1978)  Table  A.l. 

yistble  Spectrum.  Blackbody  radiance  and  hence  exi- 
tance  (emittance)  Is  negligible  in  the  visible  spectrum. 


(2.18) 


Even  at  1  fim  at  300°K  (27°C  or  80°F), 
Lilp.m,300K)  -  1.76f— 7  W/Hm2.  At  shorter 
wavelengths  (all  the  visible  wavelengths)  and/or  lower 
temperatures  (the  normal  range  of  temperature  below  the 
mesosphere)  the  emittance  is  still  less.  Hence  for  the  visi¬ 
ble  spectrum  Eq.  (2.3)  becomes 

{.•(z.0,4)  ■  Lts(z,e,d) .  (2.12) 


Equilibrium  Radiance 

The  equilibrium  radiance  is  defined  as  the  radiance 
when  the  incremental  change  in  radiance  over  the  incre¬ 
mental  path  length  dL/dr  equals  zero,  therefore  from 
Bq.  (2.1)  (Eq.  (11)  Duntley  etal.  (1957)) 

0  -  L,(z,M)  -  «(z)  I?(*,M)  (2.13) 

or  (Eq.  (40)  Chandrasekhar(1960)) 

LqizM)  -  /« (z) .  (2.14) 


Thus  the  equilibrium  radiance  is  equal  to  the  Chan¬ 
drasekhar  (1960)  source  function  J. 

No  Absorption.  In  a  medium  or  at  wavelengths  with 
no  absorption,  the  attenuation  coefficient  would  equal  the 
scattering  coefficient  and  the  equilibrium  radiance  would 
be  equal  to 

L9(s,9,4)  -  i^(i ,0,*)/z(z) .  (2.15) 

Scatttrlng  Ntgligiblt.  When  scattering  is  negligible, 
such  as  at  longer  wavelengths  in  the  infrared,  the  attenua¬ 
tion  is  equivalent  to  the  absorption.  Substituting  Eq.  (2.7) 
into  Eq.  (2.14)  we  have 

I,(*,M)-I(X,n .  (2.16) 

Equation  ofTrantftr.  An  alternate  form  of  the  equa¬ 
tion  of  transfer  is  to  divide  both  sides  by  the  attenuation 
coefficient  and  let  dr  equal  sec#  dx  [Eq.  63  Chandrasekhar 
(I960)), 

coaS  dL  (z,6,0)  /  (a(z)  dr] 

--/.(»,«,«)  + i?(*.M)  •  (2.17) 

The  cos0  is  equal  to  the  Chandrasekhar  (1960)  function  n 
and  a  Or)  dx  is  equal  to  dl(x)  the  differential  of  the  optical 
thickness  t(x). 


ISSSSnSSMJM  Length 

For  paths  of  sight  at  zenith  angles  0-0  to  70°  and 
110  to  180*,  the  incremental  path  length  is  a  simple  func¬ 
tion  of  the  altitude  increment  and  the  zenith  angle. 


Ar  -*  Ar  sec9 

When  used  for  all  paths  of  sight,  this  is  called  the  plane 
parallel  atmosphere  approximation.  The  A r  is  always 
non-negative  since  z  is  defined  as  zi-z2  (the  subscripts 
increase  with  the  flux  direction).  See  Fig.  2-1. 


Fig.  2-1.  Path  length  geometry  for  steeply  inclined  putns  of  sight. 

For  upward-looking  paths  of  sight  from  70  to  90 
degrees,  the  A r  for  a  constant  Az  is  significantly  shorter  at 
altitude  than  at  ground  level  due  to  the  curvature  of  the 
earth.  For  paths  of  sight  at  zenith  angles  from  90  to  110 
degrees,  the  A  r  for  a  constant  Az  Is  significantly  longer  at 
ground  level  than  at  altitude  also  due  to  the  curvature  of 
the  earth.  Therefore,  for  these  paths  of  sight,  the  incre¬ 
mental  path  length  A  r,  is  computed  from 

Ar,  -  [  1  -  {  nix)  ((+z)sin0  /  (»(z,)  ({+z,)l  J1  ]-^  Az  .  (2.19) 

This  is  the  classical  equation  for  computing  incremental 
path  length  at  paths  of  sight  affected  by  earth  curvature 
and  refraction.  The  nix)  is  the  refractive  index,  x  is  the 
sensor  or  observer  altitude,  £  is  the  radius  of  the  earth. 
Use  of  the  mean  earth  radius  of  6.371229E6  meters 
results  in  a  maximum  error  at  the  equator  of  -0.1  percent 
and  at  the  pole  of  +0.2  percent.  See  Fig.  2-2  for  the  rela¬ 
tionship  of  9  and  9"  for  the  downward  (and  upward)  paths 
of  sight  and  see  Duntley  etal.  (1976)  Section  2  for  a 
derivation  of  Eq.  (2.19). 


Fig.  2-2.  Path  length  geometry  for  grazing  downward  paths  of 
sight  in  refractive  spherical  atmospheres.  For  upward-looking 
paths  of  sight  the  positions  i  and  i,.  and  9  and  0  should  be 
reversed. 
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2.2  Irradiance 


Before  deriving  the  equation  of  transfer  for  irradi- 
ance,  let  us  first  define  irradiance  and  some  other  terms 
The  downwelling  irradiance  is  the  integral  of  the 

radiances  from  zenith  angles  0  to  90  degrees  weighted  by 
the  cosine  of  the  zenith  angle 

£(*,rf)  -  f  £(*,«',<*')  cosfl'  dtt  .  (2.20) 


Similarly  the  upwelling  irradiance  £(z,u)  is  computed  for 
angles  90  to  180  degrees  by  weighting  by  the  cosine  of  the 
nadir  angle  (180*0)  so  that  the  cosine  is  positive.  These 
are  the  irradlances  on  a  flat  surface  oriented  to  receive  the 
downward  and  upward  radiances  respectively. 

SflilfMlOtt 

A  second  type  of  Irradiance  is  the  scalar  or  non* 
.  directional  irradiance  in  which  the  cosine  term  is  not 
present.  The  total  scalar  irradiance  «(*)  is  defined  as 

«U)  -  f  4  li,e\dVn  .  (2.21) 


The  total  scalar  irradiance  is  equivalent  to  the  Chan* 
drasekhar  (1960)  average  Intensity  J  times  4tr.  The  scalar 
irradiance  is  related  to  the  incoming  radiant  density  w  by 

w -«<*)/  c,  (2.22) 

where  c  Is  the  speed  of  light. 


Sffltor  EftttmttJEmlUinM) 

Exitance  (omittance)  M  is  <  eflned  with  the  cosine 
term  similar  to  the  irradiance  in  Eq.  (2.20).  Let  us  define 
a  second  type  of  exitance  (omittance)  as  the  scalar  or 
non*directional  exitance  in  which  the  cosine  term  is  not 
present.  The  scalar  exitance  due  to  absorption  m(x,D  Is 
defined  as 


ei(x.n  -  L(k,T)  -4»rL(X,r)  .  (2.23) 

The  scalar  exitance  is  related  to  the  outgoing  radiant  den¬ 
sity  w(X,  T)  by 


w(k,T)  ■*  m(\,T)/c  .  (2.24) 

Let  us  also  define  another  scalar  term,  scalar  exi¬ 
tance  per  length  m*.  This  is  the  unweighted  integral  over 
4rr  of  the  path  function 

m,(t)  -  f  £,(*.«, .  (2.25) 


When  the  path  function  due  to  scaitering  is 
integrated  over  47r  it  can  be  expressed  as  a  double  integral 
by  substituting  in  Eq.  (2.4) 

m^z)  -  ^  J" /. dil  .  (2.26) 

Now  since  the  order  of  integration  is  unimportant  over  the 
double  integral,  the  relationships  in  Eq.  (2.6)  and  (2.21) 
can  be  substituted  so  that  we  have  the  scalar  exitance  per 
length  due  to  scattering  as 


-  s(z)  «(z)  .  (2.27) 

Since  the  path  function  due  to  absorption  is  isotro¬ 
pic  the  scalar  exitance  due  to  absorption  is 

m,  (z)  -  a{z)Hk,T)  4ir .  (2.28) 

O 


or 

«.  (i)-a(z)«(x,r).  (2.29) 

C 

Equation  of  Transfer  for  Irradiance 

The  equation  of  transfor  for  radiance  can  be  used  to 
obtain  an  equation  of  transfer  for  Irradiance  by  integrating 
over  4ir.  First  let  us  rewrite  Eq.  (2.1)  in  terms  of  the 
separate  components  of  absorption,  scattering,  and  emis¬ 
sion  and  express  the  incremental  path  length  in  terms  of 
Eq.  (2.18)  thus 

dL(z,0,O)cos0/<fc  »  -fl(z)£(z.M) 

-  t(i)Hz,8,<t>)  +  Lt  (z,S,<A)  +  L ,  (z)  ,  (2.30) 

$  a 

Now  multiplying  both  sides  by  </n  and  expressing  both 
sides  as  an  integral  over  4ir  we  get 


d  j  L{t,9.d>kos9dn  /di  - 

-  o(r)  J Hz,0,d>)d O  -  r(z)  |* Hi,9,<b)dft 

+  f  +  Ltg(i)  fdCl  .  (2.31) 


This  can  be  simplified  using  Eqs.  (2.20),  (2.21),  (2.27) 
and  (2.29)  to  become 

d\Eh,d)  -  £(*,«)) /eft - 

-  o(x)«(z)  -  j(z)«(r)  +  i(r)«(z)  +  o(z)m(X,D  (2.32) 


This  further  simplifies  to 

d[Eiz,d)  -  EU,u)]/dt  -  -aU)t(z)  +  o(x)  m(X,D  .  (2.33) 


secHd n  -  0  . 


(2.39) 


This  is  the  equation  for  the  net  irradiance 
E(z,d)-E(z.u)  change  with  altitude.  Let  us  define  f  (:) 
as  the  net  irradiance 


f (z)  ~  EU.d)  -  Ei2.11)  .  (2.34) 


Equation  (2.33)  can  now  be  rewritten  os 

d(UVdz -«(r)Ut(X.r) -«(:))  ,  (2,35) 

The  net  irradiance  fit)  is  equivalent  to  the  Chandrasekhar 
(1960)  net  flux  *F. 

Visible  Spectrum.  In  the  visible  spectrum  where 
Lik.T)  and  hence  mik.T)  are  negligible,  Eq.  (2.35)  can 
be  written  u  (Eq.  lO  Section  1.2  Preisendorfer  (1976) 
slnca  the  altitude  r  is  the  negative  of  the  depth  Into  the 
medium) 


d(U)  /  dz  -  -aU)*U)  .  (2.36) 

The  net  irradiance  either  stays  constant  with  altitude  Indi- 
eating  no  absorption,  or  decreases  as  depth  Into  the  atmo- 
iphere  Increases.  This  can  be  and  is  being  used  at  a  test 
for  the  Internal  consistency  of  visible  spectrum  data  when 
measured  at  a  series  of  altitudes  over  a  period  of  several 
hours. 

From  Eq.  (2.36)  we  alto  see  that  in  a  purely  scatter* 
ing  atmosphere  there  Would  be  no  change  in  net  irradiance 
with  altitude. 


Bmittoi  9f  Trtniftr  for  talir  Irrrtlinrc 

An  alternate  expression  Is  derivabla  from  Eq.  (2.30) 
if  we  flrtt  multiply  both  sidea  by  sac 9  as  well  aa  dtt  and 
than  integrate.  Thus 

rfjfi(i.a,*)rfn/4ir 

■  -fl(r)  j LU,$,6)stc$dn  ~  tU)  J L(i,t>,4)s*c9dtt 
+  fl,j(i,a,e)arcarfn  +  X,^(»)  fttcOdn  .  (2.37) 


Again  the  middle  two  terms  on  the  right  hand  side  cancel 
each  other  out  since  the  order  of  Integration  is  not  Impor¬ 
tant  in  a  double  integral 


j(i)  ^ Ut,9,4)itc9dn 


Z,(f,»\*')<r(*10,)</ni«vtldn  . 


(2.38) 


Also 


since  sec(ir-a)  -  ~str9.  Thus  the  fourth  term  also  drops 
out.  Now  Eq.  (2.37)  can  be  written  in  terms  of  the 
change  in  scalar  irradiance  with  altitude  as 

dt(:)/d2  ■  ~o(j)  J .  (2.40) 

Thus  the  change  in  scalar  irradiance  with  altitude  is  solely 
a  function  of  absorption  and  the  radiance  distribution. 
The  scattering  and  the  path  functions  due  to  scattering  and 
absorption  do  not  effect  the  change  in  scalar  irradiance  at 
any  wavelength  including  the  infrared  wavelengths. 
Therefore,  when  there  is  no  absorption  the  total  scalar 
irradiance  does  not  change  with  altitude. 

An  alternate  way  of  expressing  Eq,  (2.40)  Is  to 
separate  out  the  sun  term  and  to  express  the  more  diffuse 
radiances  as  a  radiance  difference,  sky  minus  terrain  thus 

d*U)/dz  -  -  aU)ttU)secBt 
-  aU)  f  (Z. (*,»’, tf‘)  -  LU.tr  -  8\6‘)]sec»' dl\  ,  (2,41) 


where  ,*(z)  is  the  sun  scalar  irradiance.  Equation  (2.41) 
appears  indeterminant  because  secO  goes  to  Infinity  at  90 
degrees.  However,  since  the  radiance  difference 
L(z.9)-L(z,n-9)  goes  to  zero  as  sec9  goes  to  oo,  the 
equation  Is  not  indeterminant  and  can  probably  be 
evaluated  with  measured  radiance  distributions. 

It  is  useful  at  this  point  to  speculate  on  the  relative 
importance  of  the  two  terms  In  Eq.  (2.41).  First,  we  note 
that  at  the  top  of  the  atmosphere,  the  sky  radiances  (9-0 
to  90  degrees)  are  zero  so  that  the  second  term  in  the 
equation  Is  positive  and  the  loss  due  to  absorption  at  all 
wavelengths  is  less  than  due  to  the  solar  beam  as  you 
begin  to  descend  into  the  atmosphere. 

Now,  for  additional  Insight  Into  the  two  terms  for 
the  visible  spectrum,  let  us  refer  to  our  existing  airborne 
data  on  apparent  sky  and  terrain  radiances  from  ground 
level  to  6  kilometers.  Figure  2-3  contains  two  graphs  for 
the  A  VIZ  pseudo-photoplc  filter  (see  Appendix  A  for  sen¬ 
sor  spectral  characteristics),  one  for  high  altitude,  one  for 
low  altitude,  of  measured  sky  and  terrain  radiances  for 
azimuths  from  the  sun  d>  of  0,  90,  180  and  270  degrees. 
These  graphs  are  typical  of  terrains  with  Blbedos  of 
approximately  0.1.  The  sky  radiances  (ft  -  0  to  90 
degrees)  are  generally  less  than  the  apparent  terrain  radi¬ 
ances  (9  -  90  to  180  degrees)  at  high  altitude  and  gen¬ 
erally  greater  than  the  terrain  at  low  altitude.  Thus  at  high 
altitude  the  absorption  loss  or  change  In  «(;)  with  altitude 
is  less  than  the  absorption  due  to  the  solar  beam  at  9,, 
term  1,  since  the  second  term  is  positive.  At  some  Inter¬ 
mediate  altitude  the  loss  is  equivalent  to  that  due  to  the 
solar  beam  since  the  second  term  In  Eq.  (2.41)  becomes 
zero.  At  low  altitude  the  loss  Is  greater  than  due  to  the 
solar  beam  since  the  second  term  becomes  negative 
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On  the  same  flight,  graphs  Tor  broad  band  filters  2 
and  3  with  mean  wavelengths  478  and  664  nm.  were  simi¬ 
lar  but  the  broad  band  filter  S  mean  wavelength  76S  nm 
with  an  albedo  of  0.3.  showed  sky  radiance  lower  than  ter¬ 
rain  radiance  at  low  altitude  indicating  the  absorption  loss 
would  be  leas  than  that  indicated  by  the  first  term  of 
Eq.  (2.41)  at  all  altitudes.  This  may  be  typical  of  measure¬ 
ments  for  relatively  high  albedos  and  high  transmittances. 

Normalized  functions  have  the  convenient  property 
of  always  Integrating  to  one.  Thus  from  Eq.  (2.6)  we  see 
that  the  normalized  scattering  function  is  <r<2,/3)/s(z) 
since 


[<r(z,/9)  I  s(r))  rffl  -  1  .  (2.42) 


similarly  the  normalized  scattering  path  function  from 
Eq.  (2.27)  Is  £^(*,M)/(*(f)«(2))  since 

^t,j(r,M)dn  /  (z(i)«(z))  -  I .  (2.43) 


The  single  scattering  albedo  u  is  defined  as  the  ratio  of  the 
scattering  to  the  attenuation  coefficient 

w(»)  ■  s(z)/a(z)  .  (2,46) 

The  normalized  equilibrium  radiance  can  be  expressed  in 
terms  of  the  single  scattering  albedo  as 


£f(t,M)  dn  /  |«<»)a(z)  +  (1  -oi(i)l  m(X.D)  -  1  .  (2.47) 

Vlslbk  Sptctrum,  For  the  visible  spectrum  where 
emission  is  negligible  the  integral  of  the  equilibrium  radi¬ 
ance  uivided  by  the  total  scalar  Irradiance  is  equal  to  the 
single  scattering  albedo 

(i.M)  rfn  /«(*)-  «(*) .  (2.48) 


Thus  when  there  is  no  absorption  the  equilibrium 
radiance  divided  by  the  total  scalar  irradiance  is  the  nor¬ 
malized  equilibrium  radiance, 

(r.tf.tf.)  dll  /.(D-  I  .  (2.49) 


2.3  Integrating  the  Equation  of  Transfer 

The  equation  of  transfer  for  radiance  Eq.  (2.1)  when 
integrated  over  path  length  r  yields  (Eq.  I  Duntley  et  al. 
(1957),  Eq,  1  Section  3.15  Preisendorfer  (1976)  and 
Eq.  50  Chandrasekhar  (I960)) 

-  L0 </,.», 4)  Tr(:,H)  +  ■  (2.50) 

This  equation  relates  the  apparent  radiance  L,  (2,0,0)  at 
the  sensor  altitude  2  to  the  inherent  radiance  f.„(z,,0.4)  at 
the  target  or  object  altitude  z,  to  the  radiance  transmit¬ 
tance  r,(2,0)  and  the  path  radiance  L?(z,9,<t>).  The 
transmittance  is  expressed  as  a  function  of  0  but  not  <t> 
since  transmittance  is  not  azimuth  dependent  as  long  as 
the  attenuation  coefficient  is  assumed  to  be  solely  a  func¬ 
tion  of  altitude,  the  attenuation  coefficient  is  isotropic 
horizontally. 

The  radiance  transmittance  is  defined  as 

7j,(x,0)  -exp I  1  .  (2.51) 

The  /ai(z)  is  the  optical  thickness  in  the  altitude  interval 
Az  (2,  to  z).  The  optical  thickness  Is  equal  to  the  integral 
of  the  attenuation  coefficient  from  z,  to  z 

rA,(i)  -  J*  aiz)<b  .  (2.52) 

The  path  radiance  Is  the  integral  of  the  path  func¬ 
tion  times  the  intervening  transmittance  (Duntley  el  al, 
(1957)  Eq.  (17)) 


L*(z,ti,<t>)  *’  J 9,<i>)  Tt ( (r.Wlirt'Wrfr  (2,53) 
9 

where  r,  is  from  z,  to  z. 


Horizontal  Case 

The  path  radiance  for  the  horizontal  case  can  be 
expressed  as 

90,*)  -  £.<z,90,«A)  /  Tr  1  (z,90)i/r  (2.54) 

O 

since  the  path  function  is  homogeneous  horizontally  if 
the  attenuation  coefficient  is  horizontally  isotropic.  By 


Whereas  the  normalized  path  function  for  a  medium 
which  emits  ss  well  as  scatters  according  to  Eqs.  (2.27) 
and  (2.29)  is 

f£.<x.M)rfn  /  (i(z)i(z)  +  fl(i)«.U,nl  -  I  .  (2.44) 


Since  the  equilibrium  radiance  is  the  path  function 
divided  by  the  attenuation  coefficient  (Eq.  2.14)  the  nor¬ 
malized  equilibrium  radiance  for  an  emitting  and  scatter¬ 
ing  medium  is 


L, (z.Mlrfft  a(z)/|s(z)«(*)  +  a(z)mti,T)]  -  I  .  (2.45) 
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multiplying  the  numerator  and  denominator  by  aiz)  we 
cen  integrate  to  get 

£,*<j,90, *>  -  £.(x.90,4)(l  -  r,(z, 90))  /  a(z)  (2,55) 


or  from  Eq.  (2.14) 


£,*<», 90.4)  -  L,(t,9 0.4)  II -r,(z.90)J  .  (2.56) 

This  i»  the  atandard  equation  for  the  horizontal  path  of 
eight  which  holds  for  a  scattering,  absorbing  and  emitting 
medium  which  ia  horizontally  homogeneous. 

Since  the  sky  radiance  is  the  path  radiance  from 
space-to-eirth,  the  horizon  sky  radiance  is 


Ll  (|,90,4)  -  l, <1,90,4)  ll  -  r.(i,90)l  .  (2.57) 

The  transmittance  from  space-to-earth  at  zenith  angle  9  is 
a  function  of  the  vertical  transmittance  and  the  relative 
optical  airmass  mM(r,0)/m.(z,0).  Therefore 

Tjt, 90)  -  r«(i,0)"-(,'W)/"“(,  0> .  (2.58) 


Air  Maas 

The  absolute  optical  alrmasa  mm(z,9)  is  the  mass  of 
an  aircolumn  of  unit  cross  section  which  is  the  integral  of 
the  density  over  the  path  length  from  altitude  i  to  «> 


mM(i,e)  -  p(t)dr  .  (2.59) 


The  relative  optical  airmass  is  the  ratio  of  the  airmass  at 
some  angle  9  to  the  vertical  airmass  m.(z,0)/m«(z,O),  In 
computations  of  the  optical  air  mass,  the  incremental  path 
length  A r,  is  computed  using  Eq.  (2.19)  which  takes  into 
account  both  the  curvature  of  the  earth  and  refraction. 
For  the  computations  herein  the  mean  earth  radius  was 
used. 

The  square  of  the  refractive  Index  is  expressed  as  a 
function  of  the  density  ratio  by  Kasten  (1965) 

(*<*)/ s(ff))1  -  1  +  2|*(r)  -  l]|l  -  pitt)  /  pit)]  (2.60) 


To  a  good  approximation,  the  refractive  index  can  be 
expressed  as  a  function  of  density  as 

|«(t)  -  Ij  /  pit)  -  l*<0)  -  l)  /  p<0)  ,  (2.61) 

Substituting  Eq.  (2.61)  into  Eq.  (2.60)  we  get  an  equation 
in  terms  of  the  refractive  Index  at  sea  level  as  follows: 


ln(z)/  ni:^)1 

-  1+2  [rt(0)  -  l)|p(z)//>(0)  -  p(:,)/p( 0)]  .  (2.62) 

The  refractive  index  is  determined  using  the  disper¬ 
sion  formula  of  Edlen  (1953)  for  the  refractive  modulus 
/V(0)  at  15  degrees  Celsius  sea  level  appropriate  for  a  den¬ 
sity  of  1.2250  kg  m~i  (the  standard  atmosphere  sea  level 
density), 

M0)  -  64.328  +  29498.10/  (146  —  1/\J) 

+  255.40/(41-  I  /*»).  (2.63) 

where  k  is  in  pm  The  refractive  modulus  is  related  to 
the  Index  of  refraction  by 

/v(0)-lrt(0)-l)l£6.  (2.64) 

Computations  of  relative  airmass  were  made  for  sea 
level  and  6  kilometers  altitude  using  density  values  based 
on  the  V.  S.  Standard  Atmosphere  (1976)  and  the  warmest 
and  coldest  supplementary  atmospheres  in  the  U.  S.  Stan • 
dard  Atmosphtrt  Supplements  (1966).  The  density  at  each 
altitude  was  obtained  by  truncated  Chebyshev  expansion 
using  coefficients  for  the  atmospheres  between  0  and  80 
kilometers  U.  S.  Standard  Atmosphere  Supplements  (1966) 
Tables  4.2,  4.3,  4.6,  and  4.11).  Computations  for  the 
stsndard  atmosphere  tor  0  to  200  kilometers  agreed  within 
1  percent  with  the  results  for  0  to  80  kilometers,  therefore 
the  results  reported  in  Table  2.1  are  for  0  to  80  kilometers 
for  all  three  atmospheres. 

The  computations  were  made  for  a  wavelength  of 
700  nanometers.  Oast  et  al.  (1965)  indicate  less  than  a  1 
percent  change  in  relative  air  mass  for  even  a  90  degree 
zenith  angle  for  wavelengths  from  300  to  3000  nanome¬ 
ters,  The  values  in  Table  2.1  for  the  standard  atmosphere 
at  sea  level  agree  within  0.1  percent  of  the  Oast  etal. 
(1965)  values  for  700  nanometers  except  at  90  degrees 
where  they  agree  within  0.9  percent. 

Table  2.1  presents  values  of  relative  optical  airmass 
for  angles  70  to  90  degrees.  From  0  to  70  degrees,  the 
relative  airmass  is  equal  to  sec9  within  1  percent.  The  45 
degree  North  Spring/  Fall  values  are  appropriate  for  the 
U.  S.  Standard  Atmosphere  (1976).  For  purposes  of  com¬ 
parison,  the  values  for  Bemporad  from  Cast  etal.  (1965) 
are  given  in  Column  2  and  for  Kasten  (1965)  in  Column 
3.  Bemporad  assumed  a  sea  level  temperature  of  0 
degrees  Celsius  and  a  wavelength  of  540  nanometers. 
Kasten  used  the  ARDC  (1959)  densities  which  are  essen¬ 
tially  equivalent  to  the  U.  S.  Standard  Atmosphere  (1976) 
and  a  wavelength  of  700  nanometers. 

Our  computations  compare  well  to  the  Kasten 
(1965)  values  except  at  zenith  angle  90  degrees.  Since  his 
90  degree  value  is  also  at  odds  with  the  Oast  etal.  (1965) 
values,  his  90  degree  approximation  technique  may  be 
inadequate  (there  are  special  computational  problems  at  90 
degrees).  Bemporad’s  values  lie  as  expected  between  the 
45  degree  Spring/Fall  values  and  the  60  degree  North 
(cold)  values  as  does  his  ground  level  temperature. 


Table  2.1.  Relative  optical  uir  mass. 


In  ordsr  to  evaluate  tho  effect  of  aenior  altitude  on 
th«  relative  air  maas,  the  6  kilometer  values  are  divided  by 
the  sea  level  values  and  the  ratio  given  in  Table  2.2  for  86 
to  90  degrees.  The  relative  airmass  values  agree  for  6 
kilometers  with  the  sea  level  values  within  1  percent  from 
70  to  86  degrees  regardless  of  temperature.  The  largest 
altitude  difference  is  +4.1  percent  for  the  US.  Standard 
Atmosphore  at  90  degrees  zenith  angle.  Altitude  affects 
the  relative  airmass  less  than  1  percent  for  the  coldest 


atmosphere  except  at  90  degrees  (the  temperature 
difference  0  to  6  kilometers  is  also  rmallest  for  the  60 
degree  North  (cold)  atmosphere). 

The  effect  of  temperature  on  the  relative  airmass  is 
evaluated  by  ratioing  the  warmer  and  colder  atmosphere 
values  to  the  Standard  Atmosphere  values  at  45  degrees 
North  in  Table  2.2.  The  effect  of  temperature  Is  generally 
greater  than  the  effect  of  sensor  altitude.  Also  the  effect 
of  temperature  Is  greater  at  sea  level  than  at  6  kilometers. 


Table  2.2.  Altitude  end  temperature  effects  upon  relative  optical  airmass. 


Horizon  Sky  Equilibrium  Radiance 

As  can  be  seen  from  Eq.  (2.57)  when  the  transmit¬ 
tance  r.(z.90)  is  negligible,  the  horizon  sky  radiance  is 
equal  to  the  equilibrium  radiance.  In  order  to  determine 
the  minimum  vertical  transmittance  necessary  to  produce 
a  horizon  sky  radiance  at  equilibrium  let  us  rearrange 
Eq.  (2.58)  as  follows 

In  rjt.O)  -  In  r„U.%)/[«o.<i.90)/m..(zl0)l .  (2.65) 


For  3  percent  accuracy,  f.(z,90)  would  be  0.03,  then 
solving  Eq.  (2.(5)  using  the  relative  airmass  value  for  the 
standard  atmosphere,  we  And  that  7*M(z,0)  -  0,911.  Thus 
for  all  vertical  transmlttances  0.911  or  less,  the  horizon 
sky  will  be  the  equilibrium  radiance  to  an  accuracy  of  3 
percent  or  better.  The  Rayleigh  transmittance  for  a 
wavelength  of  555  nanometers  for  (he  Standard  Atmo¬ 
sphere  Is  0.910.  Thus  for  all  wavelengths  555  nanometers 
or  smaller,  the  horizon  sky  radiance  will  always  be  equal 
to  the  equilibrium  radiance  at  sea  level.  Since  the  photo- 
pic  haa  been  found  to  be  similar  to  the  narrow  band  meas¬ 
urements  at  555  nanometer,  ibe  photoplc  horizon  sky 
radiance  at  sea  level  is  also  «lw«y;v  reasonably  equal  to  the 
equilibrium  radiance  for  non  days. 

The  Rayleigh  horizontal  transmittance  at  6  kilome¬ 
ters  can  be  related  to  the  vertical  sea  level  transmittance 
by 

*  r„(6,90)  -  ft  rM<0,0)"-(‘  ’OI/"-(0,01  .  (2.66) 

Values  of  absolute  vertical  alrmaaa,  and  relative  vertical 
airmass  for  sea  level  and  6  kilometers  are  given  in  Table 
2.3  for  the  same  three  atmospheres  as  in  Table  2.1.  Since 


percent.  This  turns  out  to  be  r.„( 0,0)  -  0.826.  Since  the 
Rayleigh  transmittance  for  465  nanometers  wavelength  for 
the  Standard  Atmosphere  is  0.823,  the  horizon  sky  radi¬ 
ances  at  6  kilometers  for  wavelengths  465  nanometers  or 
less  will  always  be  at  equilibrium  in  cloudless  skies. 

2.4  Integrating  the  Equation  of  Transfer 
for  Irradiance 

The  equation  of  transfer  for  irradiance  Eq.  (2.35) 
can  be  rearranged  so  that  both  sides  can  be  integrated  with 
respect  to  altitude.  The  limits  of  Integration  are  between 
z„  and  t  as  we  descend  into  the  atmosphere  so  that  z„  is 
greater  than  2, 

j  dfit)  -  f  a{i)[m(\,T)  -  .  (2.68) 

The  left  hand  term  is  easily  integrated  but  the  right  hand 
term  is  not.  Partially  integrated  Eq.  (2.68)  becomes 

f(t)  -- f(r0)  - JafzriMta.n -«(i)U.  (2.69) 


Approximate  Solution 

A  reasonable  approximation  to  the  integration  can 
be  obtained  If  we  assume  that  the  change  in  the  difference 
t  m(k,T)-*{t) )  is  small  enough  from  z0  to  x  such  thst 
»n  average  value  { m(k,  )  can  be  substituted. 

Substituting  in  the  average  scalar  Irradiance  and  exltance 
difference  which  can  now  be  taken  out  of  the  integral  we 
get 


m„(6,90)/«.(0.0) 


m,(6,90)  w,(6,0) 

id,  (6.0)  mJojST  ' 


(2.67) 


(U)  -  f(*0)  -  I  n(k, T)  - 


a(t)dt  . 


(2.70) 


We  now  have  the  bails  for  evaluating  Eq,  (2.66). 
Using  an  equation  similar  to  Eq.  (2.65),  we  can  now 
evaluate  the  vertical  Rayleigh  sea  level  transmittance 
necessary  to  produce  a  horizon  sky  radiance  at  6  kilome¬ 
ters  squal  to  the  equilibrium  radiance  to  an  accuracy  of  3 


The  integral  of  the  absorption  with  altitude  is  the  optical 
thickness  due  to  absorption  ,  rA,  (r)  therefore 

fU>  -  f(j„)  -  #r4,(z)  I  «?a77V~' T7S5T 1 .  (2.71) 


Table  2.3.  Temporal  and  latitudinal  comparisons  of  vertical  optical  airmass. 
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The  error  in  the  resultant  optical  thickness  is  probably  less 
than  the  variability  or  the  difference  lm(A,  D-«(r)]  from 
the  average  value.  This  is  an  approximate  equation  for 
the  net  irradiance  change  with  altitude. 

An  expression  for  the  absorption  optical  thickness 

can  be  obtained  by  rearranging  Eq.  (2.71) 

a'u  -  I  1  >  I  "tTTTr-  atari  )  .  (2.72) 


Thus  if  one  has  a  measure  of  ambient  temperature  and 
downwelling,  upwelling  and  scalar  irradiance  at  two  alti¬ 
tudes  it  and  t  where  z,  is  greater  than  z,  the  absorption 
optical  thickness  can  be  obtained. 

No  Absorption,  When  there  is  no  absorption  and 
hence  no  emission,  Eq.  (2.71)  becomes 


{<*)  0  (2.73) 

or  substituting  in  Eq.  (2.34) 

£(z,rf)  -  £■(*,«)  -  £(z„,d)  -  £<z0,u) .  (2.74) 


Whan  there  is  no  absorption  the  net  irradiance  does  not 
change  with  altitude. 

The  albedo  A(t)  is  defined  as  the  ratio  of  the 
upweillng  to  the  downwelling  irradiance. 


A(s)-£(t,u) /£(»,</).  (2,75) 


Substituting  Eq.  (2.73)  into  Eq.  (2.74)  we  get  an  expres¬ 
sion  in  terms  of  downwelling  irradiance  and  the  albedo 

£<z,rf)|l  -  A  (i)l  -  £<z0,rf)(  1  -  A  (*„))  .  (2,76) 


Thus  the  apperent  albedo  change  with  altitude  is  directly 
linked  to  the  downwelling  irradiance  change  with  altitude 
when  there  is  no  absorption. 

The  albedo  out-of-the-atmosphere  A  (*>)  can  be 
computed  when  the  albedo  and  the  downwelling  irradiance 
at  altitude  i  are  known  since  the  downwelling  irradiance 
out  of  the  atmosphere  la  the  solar  scalar  irradiance  ,«(<») 
times  the  cosine  of  the  sun  zenith  angle 

£(«M)  -  ^Mcoitj  .  (2,77) 


Substituting  Eq.  (2.77)  into  Eq.  (2.76)  and  rearranging  we 

A(~)  -  1  -  £(z,rf)|l  -  /4(z)l  /  l,«McoiS,)  .  (2,78) 

For  example,  if  we  assume  the  photoplc  clear  day  to 
have  negligible  absorption,  we  can  use  the  illuminances 


from  Brown  (1952)  to  compute  the  albedo  out  of  the 
atmosphere  assuming  various  albedos  at  ground  level. 
Using  ,«(«')  -  1.89/T3  W/m}nm  and  the  Brown  (1952) 
illuminances  converted  to  irradiunce  in  the  same  units 
(see  Appendix  A  for  details  on  the  photopic  as  a  broad¬ 
band  sensor),  the  results  urc  shown  in  Table  2.4  and  Fig 
2-4  for  various  ground  level  albedos.  Ground  level  albedo 
A  (0)  equal  to  0. 1  is  probably  the  most  reasonable  in  rela¬ 
tion  to  the  Brown  (1952)  illuminances. 


Table  2.4.  Albedo  out-of-thc-atmospherc  for  u  non-absorbing 
atmosphere  based  on  Brown  (1952)  and  Eq.  (2.78). 
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Fl«.  2-4.  Albedo  out-of-the-aimospherc  for  a  non-absorbing 
atmosphere,  based  upon  Brown  (1952)  and  Eq. (2.78). 


Integrating  for  Scalar  Irradiance 

The  equation  of  transfer  for  scalar  irradiance 
Eq,  2.40  can  be  Integrated  by  first  dividing  both  sides  by 
«(z).  Now  multiplying  by  dz  and  integrating  from  z„  to  z 
where  za  is  greater  than  z  as  we  descend  into  the  atmo¬ 
sphere, 


</«(z)/«(z) 


«<z) 


snUdildi 


(2.79) 


-II- 


Again  the  left  hand  term  is  easily  integrated  but  the  right 
hand  is  not.  Let  us  define  a  function  D(z)  which  is  a 
semi-normalized  function  of  the  4t r  radiance  distribution 
at  altitude  z 


H:,8,<j>)sec9dl\  /  <(2)  . 


The  D(z)  is  the  unitless  weighting  function  for  the 
integral  of  absorption  with  altitude  as  it  affects  the  change 
of  scalar  irradiance  with  altitude.  An  alternate  way  of 
expressing  Eq.  (2.80)  is  to  separate  out  the  sun  term  and 
to  express  the  more  dilftise  radiances  as  a  radiance 
difference,  sky  minus  terrain  as  in  Eq.  (2.41) 

0(f)  -  step,  t*{z)/*{z) 

+  J" -  L(:,n  -  seed'd  n  /  «(r)  .(2,81) 

The  first  term  is  exact  but  the  second  term  at  large 
zenith  angles  becomes  an  approximation  since  stcO' 
becomes  very  large  near  90  degrees.  Fortunately,  at  90 
degrees,  AL  Is  zero. 

Now  substituting  Eq.  (2.80)  Into  Eq,  (2,79)  and 
integrating  the  lefthand  side  we  get 

In  [<(s)/«(i9)]  -  -J" a(i)0(i)A  .  (2.82) 


Equation  (2.85)  indicates  that  one  needs  a  measure¬ 
ment  of  the  4rr  radiance  distribution  at  both  altitudes  in 
order  to  compute  the  scalar  irradiance  and  the  radiance 
distribution  function  D(z). 

Water  as  a  Scattering  and  Absorbing  Medium.  As  an 
interesting  and  illustrative  aside,  one  may  consider  the 
application  of  the  foregoing  equations  within  the  hydrolog¬ 
ical  context.  Although  the  principal  emission  mechanism 
in  water  is  thermal  there  is  also  bioluminescence. 
Bioluminescence  is  essentially  chemiluminescence,  that  Is, 
energy  from  chemical  reaction  instituted  by  biological 
organisms.  It  Is  also  isotropic.  Thus  in  water  the  total 
emittance  path  function  is  L»t(.d)~L»a(d)+L,b(d),  In 
water  the  radiances  are  denoted  in  terms  of  the  depth  d 
beneath  the  water  surface  rather  than  the  height  above  the 
sea  floor  z.  Since  the  path  function  due  to  biolumines¬ 
cence  is  isotropic,  the  equation  of  transfer  for  scalar  irradi¬ 
ance  Eq.  (2.40)  is  appropriate  for  water  whether 
bioluminescence  Is  present  or  not,  On  the  other  hand,  the 
equation  of  transfer  for  irradiance  Eq.  (2.35)  only  applies 
to  water  when  bioluminescence  is  not  present,  otherwise, 
a  bloluminescent  term  must  be  added.  Also  in  the  visible 
spectrum  in  water,  thermal  emission  Is  negligible  as  it  is  in 
the  troposphere.  In  water,  below  a  certain  depth  d\  the 
light  distribution  becomes  asymptotic,  with  the  sun  no 
longer  a  separate  component.  Thus  the  function  D{d\) 
which  Is  solely  a  function  of  the  lighting  distribution 
becomes  a  constant  and  can  be  expressed  as  the  right 
hand  term  of  Eq.  (2,81) 


This  ii  the  partially  integrated  equation  for  the  change  of 
scalar  irradiance  with  altitude. 

Approximatt  Solution ,  A  reasonable  approximation 
to  the  Integration  can  be  obtained  if  we  assume  that  the 
change  In  D(t)  is  small  between  z,  and  t  ao  that  an  aver¬ 
age  value  J5(Az)  can  be  substituted 


-  fttWJ)  -  L(rf„ir  -  e\*')]»c9’rfn  /  «(</,).  (2.86) 

Also  at  depth  d\  and  below,  absorption  Is  constant. 
Therefore  the  equation  of  transfer  for  scalar  irradiance 
Eq.  (2.79)  becomes 


(<(t)/«(z0)]  -  -B(Az)  fahiib. 


dt(d)/t(d)  -  -o  Did,)  I  dd  .  (2.87) 


or  in  terms  of  the  absorption  optical  thickness 

In  l«(f)/«(fe)l  -  -B(Ai)  flr4l(f)  ■  (2,84) 

The  error  in  the  resultant  optical  thickness  is  again  prob¬ 
ably  less  than  the  variability  of  the  D(z)  from  the  average. 
This  Is  an  approximate  solution  for  the  change  of  scalar 
irradiance  with  altitude. 

In  order  to  obtain  a  measure  of  absorption  optical 
thickness  between  two  altitudes,  we  would  rearrange 
Bq.  (2.84)  thus 

.^(i)  -  In  U<*0) /«(i)l/ 5(4*) .  (2.85) 


This  is  easily  integrated  with  depth  to  become 

In  [«((y)/«(<y0»  -  -  0  D(d,)(d  -  d0) ,  (2.88) 

where  d,  must  be  at  or  below  d\  the  depth  at  which  the 
lighting  distribution  becomes  asymptotic. 

In  water,  the  constant  k  is  defined  as  [Eq.  7  Section 
1,2  Vol.  I  Prelsendorfer  (1976)1 


lnM<r)/«(</p))  -kld-d0) 


Therefore  Eq.  (2.88)  together  with  Eq.  (2.89)  give  us  an 
equation  for  k  in  terms  of  the  absorption  a  and  the 


-V.  v 


-  ■>>_v  .v, 


asymptotic  lighting  distribution  function  D  defined  by 

Eq.  (2.86). 

k  -  aD(d\)  (2  90) 

Thus,  having  noted  the  conceptual  similarities  between  the 

hydrosphere  and  atmosphere,  our  discussion  returns  to 
the  subject  at  hand,  the  atmosphere  as  a  radiative  transfer 

medium. 


Discussion 

Equations  (2.71)  and  (2.84)  which  describe  net  irra- 
diance  and  scalar  Irradlance  are  essentially  three  dimen¬ 
sional.  They  address  the  effect  of  absorption  over  an  area 
as  well  as  vertically.  They  do  not  assume  homogeneity 
but  they  do  assume  a  relatively  small  change  in 
m(ksT)-t(z ),  and  D(z)  over  the  altitude  interval.  Thus 
clouds  may  be  present  and  the  equations  can  still  be  valid. 
The  equations  assume  ze  is  vertically  above  z  and  describe 
a  single  point  in  time. 

When  there  are  clouds  in  the  altitude  Increment  dr, 
Eq.  (2.85)  should  indicate  absorption.  This  Is  an  absorp¬ 
tion  based  on  the  4n  lighting  distribution.  Thus,  even  if 
there  is  no  absorption  in  the  direction  of  the  sun,  the 
presence  of  clouds  will  affect  the  4ir  lighting  distribution 
and  absorption  should  be  indicated. 

We  now  have  two  equations  Eq.  (2.72)  and  (2.84) 
to  obtain  the  absorption  optical  thickness.  Therefore  it 
should  be  possible  to  use  them  simultaneously  us  a  check 
on  the  efficacy  of  the  approximate  solutions,  that  is 

fl/4l(r)  -  I  t(z)  -  f(r„)  1  /  I  m(\,T)  -  «(Az)  1 

-lnl«(r0)/.(r)]/B(dr)  .  (2.91) 

Visible  Spectrum.  In  the  visible  spectrum,  the  emis¬ 
sion  is  negligible  therefore  Eq.  (2.91)  simplifies  to 

„'4,(*>-  l<(*0M(i)]/i<dr) 

-In  U(z0)/«(z)]/ZS<Az> .  (2.92) 

The  above  equations  are  monochromatic.  They  are, 
however,  reasonable  engineering  approximations  for  a 
broadband  »vnsor  when  the  absorption  is  primarily  In 
broadband  continua  or  due  to  the  suspended  particles  or 
water  droplets  and  not  line  or  band  spectra.  This  holds 
reasonably  woll  in  the  atmosphere  In  the  visible  portion  of 
the  spectrum. 

To  test  Eq.  (2.92)  and  to  get  an  understanding  of 
the  lighting  distribution  function  D(z),  we  used  the  sky 
and  terrain  radiances  measured  for  Filter  2  mean 
wavelength  478  nanometers  on  flight  C-466  in  Moppen, 
Germany  on  15  August  1978  during  project  OPAQUE  V 
described  in  Johnson  and  Gordon  (1980).  The  altitude  r„ 
was  chosen  as  the  highest  altitude  flown  6  kilometers  and 
the  z  used  was  for  the  lowest  altitude  flown  0.2  kilome¬ 
ters.  Thus  Az  equals  S.8  kilometers.  The  sun  was  unob¬ 


scured  at  both  altitudes  but  there  were  scattered  cirrus 

clouds. 

The  first  step  was  to  correct  the  scanner  data  in  the 
region  0  to  25  degrees  scattering  angle  from  the  sun. 
Recent  scanner  tests  indicate  spurious  sun  reflection  in  sky 
radiance  measurements  near  the  sun.  The  method  used 
was  similar  to  the  Barteneva  (I960)  method  for  extrapo¬ 
lating  the  volume  scattering  function  near  scattering  angle 
0  degrees.  We  assumed  log  sky  radiance  at  each  zenith 
angle  to  be  linear  with  cos/3  from  0  to  30  degrees.  The 
slope  of  the  line  was  established  by  the  values  near  30 
degrees  scattering  angle. 

The  next  step  was  to  obtain  a  measure  of  radiance 
transmittance  appropriate  for  the  unobscured  sun.  This 
was  done  using  the  solar  almucantar  sky  radiances  in  a 
method  described  later  in  Section  4.6. 

As  a  check,  the  resultant  computed  downwelling 
irradiances  were  compared  to  the  irradiometer  measure¬ 
ments.  The  irradiances  computed  from  the  sky-terrain 
scanner  measurements  are  given  in  Table  2.5  together 
with  the  ratios  to  the  irradiometer  values.  The  computed 
downwelling  irradiances  compare  as  well  or  belter  than 
the  upwelllng  Irradiances  and  thus  the  corrected  sky  radi¬ 
ances  were  considered  reasonable. 

For  computing  the  radiance  distribution  function 
Diz),  Eq.  (2.81)  was  rewritten  as 

D(z)  -  src9s  sii,z)/t(z) 

+  2rr  /  [Liz, 9')  -  Liz, it  -  0')]wc8W  /  «(r)  .  (2.93) 

O 

The  average  radiance  1(2.9')  for  each  zenith  angle  was 
computed  by  summing  the  radiances  for  all  60  azimuths 
and  dividing  by  60.  The  resultant  average  radiance 
differences  L(:,9')-I(z,ir-6')  are  graphed  In  Fig.  2-5. 

The  computations  of  the  distribution  function  D(z) 
using  Eq.  (2.93)  are  given  in  Table  2.6.  Note  that  the 
second  or  sky-terrain  term  is  relatively  small  (23  percent 
of  the  total)  at  high  altitude  but  more  than  half  the  total  at 
low  altitude. 

The  evaluation  of  Eq.  (2.92)  is  also  shown  in 
Table  2.6  together  with  the  average  absorption  optical 
thickness  based  on  an  average  of  the  values  using  the  two 
methods.  The  average  differs  from  the  individual  evalua¬ 
tions  by  ±9  percent.  This  compares  well  to  the  estimate 
of  error  for  each  value.  The  estimated  error  for 
Eq.  (2.72)  is  ±8  percent  based  on  the  variability  of  the 
difference  lm(k,7')-«(r))  from  the  average.  Whereas  the 
estimated  error  for  Eq.  (2.85)  is  14.7%  based  on  the  varia¬ 
bility  of  0(Az)  from  the  average  value. 

The  resultant  absorption  optical  thickness  probably 
errs  on  the  high  side  since  no  attempt  was  made  to  correct 
for  the  small  sun  zenith  angle  change  between  the  meas¬ 
urements  at  6  kilometers  and  0.2  kilometers. 

The  above  is  too  small  a  sample  to  validate  the 
method  but  is  encouraging  enough  to  indicate  further 
development  might  prove  fruitful. 
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Table  2.5a.  Compuleti  scalur  irradianees  lor  flight  C-406, 
filler  2  A  —  478/i/m. 


COMPUTED  SCALAR  IRRADIANCI  S  I 

SPECIFIC  ]  Avcrjjjc  ,  Chunitc 


IRRADIANCE 

COMPONENT 

Al  6km 

f(6) 

Al  2krr 

«(2) 

Over 

AUltucJe 

•  (A*) 

Over 

Altitude  I 

Ink  {/0i/i  (;)1  1 

Sun  ,«<:) 

1796 

1113 

| 

Skyrfe<.-.d) 

131 

its 

i 

Upwolllnj  « (j.u) 

«4 

237 

Tout  «(»> 

2331 

2171 

2361 

1601 

Alberto  ,A 

.199 

119 

1 

Table  2.5b.  Computed  flat  plate  Irradianees  for  flight  C-466. 
filter  2,  K  -  478  nnt. 


SPECIFIC 

irradiance  component 

COMPUTED  Ft.AT  PLATE  IRRADIANCES 

Al  6  km  £(6) 

Calculated  from  Scanner 

Ai  2km  ft.2) 

Calculated  from  Scanner 

Sun  ,£(*> 

1339 

737 

Sky  de<t.<n 

114 

310 

Tout  Downwellina  £(1.41 

14)3 

1137 

Tout  Upwelllrt|  Etl.k) 

117 

77 

Albedo.  All) 

127 

069 

Net  <01 

1216 

1060 

Calculaled-Maaeured  Rilloe 

Downwelllni  (Scan/lrred) 

1.04 

1.02 

Upwelllnj  (Stan/lrrid) 

096 

III 

Table  2.6a.  Lighting  dlitrlbutlon  function  D(z)  from  Eq.  (2.93), 
flight  C-466,  filter  2,  X  -  478  m, 


Table  2.6b.  Comparison  of  absorption 
optlcul  thickness  determinations. 


Sped  Ao 

DISTRIIUTION  FUNCTION  D(t) 

Dtalrlbuilon 

Al  6km 

Ai  .2  km 

Component 

Dtt) 

D(2) 

Sun 

0  921 

0.731 

Sky-Terrain 

0.274 

0163 

Total 

1.202 

1616 

Averaaa  0(4e> 

1.409 

Computational 

Procedure 

Compuied  Abeorpnon 
Optical  Thlckiiaet 

p's  »<•» 

From  Net  Irradiance 

Eq.  (2.72) 

0  0954 

Prom  Sc*l»f  Irradiance 
Eq.  (2  8$) 

0  1141 

Avereie  Value 

0.1048 

«  A 


I1|.  2-S.  Average  radiance  difference  Z(i,6)-I(i,ir-8)  for 
(light  C-466,  filler  2  (X  -  478  nm). 


Earth  Curviture  Effects 

The  above  equations  were  developed  for  a  plane 
parallel  atmosphere.  Refraction  and  earth  curvature 
effects  were  Ignored.  However,  for  upward  and  downward 
paths  of  sight  between  6  kilometers  and  ground  level,  the 
incremental  path  length  Ar,  is  affected  to  a  significant 
extent  by  earth  curvature  and  refraction  for  the  longer 
slant  paths  of  sight.  For  these  paths  of  sight,  the  incre¬ 
mental  path  length  Ar  is  computed  from  Eq.  (2.19). 

As  indicated  earlier,  the  A r  for  a  constant  A;  is 
shorter  at  altitude  than  at  ground  level  for  upward  paths  at 
zenith  angles  greater  than  70  degrees,  and  longer  at 
ground  level  than  at  altitude  for  downward  paths  of  sight 
less  than  110  degrees  in  zenith  angle.  The  ratio  of  the 
Incremental  path  length  at  ground  level  and  6  kilometers 
for  these  paths  is  given  in  Table  2.7.  The  A r  values  were 
computed  from  Eq.  (2.19)  for  a  refractive  Index  appropri¬ 
ate  for  a  wavelength  of  700  nm  and  standard  atmosphere 
densities.  There  is  no  ratio  for  zenith  angles  92  and  91 
degrees  since  at  these  angles  the  path  of  sight  does  not 
Intersect  the  surface  of  the  earth  due  to  earth  curvature. 

The  A r  at  the  sensor  position  is  equivalent  to  secO. 
It  is  the  A r  at  the  beginning  of  the  flux  path  that  departs 
from  mc0  by  the  ratios  indicated  in  Table  2.7.  These 
ratios  indicate  the  maximum  error  In  using  s ec9  for  these 


Table  2.7.  Ratio  of  incremental  path  length  Ar(:,)/Ar(c)  for 
-  30  meter  increments  for  both  upward  and  downward  paths 
of  sight  between  ground  level  (/  **0)  and  6  km  (/  -  61 


!  I 

/xrnitli  •  Upward  l’jlh*  '  Zenith  i  Downward  Pulhs  I 

Angle  cl  Sight  Angle  of  Sight 


H 

(degree*) 

■1M6I/AM0) 

'' 

(degree*! 

Ar(0)/Ar(6) 

70 

.994 

no 

1.006 

73 

919 

105 

1.012 

76 

917 

104 

1.013 

77 

913 

103 

1.016 

71 

912 

102 

1.019 

79 

979 

101 

1.022 

10 

973 

100 

1  027 

II 

969 

99 

1.034 

11 

.961 

91 

1  044 

1) 

.930 

97 

1.039 

14 

.934 

96 

1.013 

13 

901 

95 

1.127 

16 

166 

94 

1  225 

17 

792 

93 

1  369 

II 

.654 

19 

397 

angles  for  paths  between  0  and  6  km.  The  error  is  less 
thBn  3  percent  for  paths  from  70  to  80  and  100  to  110 
degrees. 

For  paths  from  80  to  100  degrees,  It  is  appropriate 
to  compute  the  relative  alrmass  to  compare  to  sec6.  The 
absolute  airmass  was  computed  using  Eq.  (2. 59)  integrat¬ 
ing  from  0  to  6  km.  Values  of  relative  airmass,  and  the 
ratios  to  sec9  for  selected  zenith  angles  from  80  to  100 
degrees  are  given  in  Table  2.8.  Use  of  seed  to  obtain  A r 
for  a  6  kilometer  altitude  Increment  introduces  an  error  of 
+  1.1  to  +4,5  percent  for  80  to  85  degrees  zenith  angle 
and  an  error  of  -1.5  to  -6.3  percent  for  95  to  100  degrees. 
This  is  small  enough  to  be  neglected.  However,  at  87 
degrees  the  error  Is  +11.6  percent  and  at  93  degrees  the 
error  is  -20,2  percent. 

Checking  back  Into  the  computations  we  found  that 
the  contribution  to  D(z)  at  zenith  angle  87.5  degrees  is 
large  at  high  altitude,  essentially  all  of  the  second  term. 
Whereas  at  low  altitude  the  radiance  difference  at  87.5 
degrees  is  less  relative  to  the  other  zenith  angles  so  that  it 
comprises  only  18  percent  of  the  second  term  or  10  per¬ 
cent  of  the  total.  The  contribution  of  the  radiance 
difference  at  87.5  degrees  Is  probably  the  most  question¬ 
able  part  of  the  evaluation  of  the  distribution  function 
D(i).  Thus  the  uncertainty  of  Dibkm)  is  on  the  order  of 
23  percent  and  of  D(0.2k.m)  is  10  percent. 
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Table  2.8.  Comparison  of  relative  optical  airmass,  mr{2,0)/mr(zt 0) 
to  seci)  (computed  using  &2  -  30  meters  for  vertical  path  of  6  km). 


upward  paths  of  sight 


Zennh 

An|le 

«, 

(deireeil 

SK» 

Relative 
Alrmtu 
m,(0 .«> 

Ratio 

_ wc» 

Rtl.  Airmass 

Zenith 

Angle 

", 

(degrees) 

<0 

5.759 

5.494 

1.011 

IX 

IS 

11.474 

10913 

1045 

95 

17 

19.107 

17.121 

1.116 

93 

downward  paths  of  sight 


Relative 

Airmen  Ratio 


mj(6, 1 80)  Rtl.  A i 'man 


3.  SOLAR  ALMUCANTAR 
3.1  Sky  Radiance 

The  almucantar  is  the  part  of  the  sky  at  a  constant 
zenith  angle,  all  azimuths.  The  solar  almucantar  is  the  sky 
at  the  zenith  angle  of  the  sun  6,.  The  sky  radiance  in  the 
solar  almucantar  is  the  path  radiance  out-of-the- 
atmosphere  to  the  sensor  at  altitude  z.  The  sky  radiance 
can  be  expressed  as 


The  sun  scalar  irradiance  and  the  transmittance  can 
be  expressed  as  a  function  of  the  sun  irradiance  out  of  the 
atmosphere  ,«(<»)  as  follows 

st(zl)Trt{z,9s)  -  f «(«o)  exp  ( —  J  a(z)dr~  J*  a(z)rfr)  .  (3.6) 


The  two  integrals  can  be  combined  and  expressed  in  terms 
of  the  limits  »  to  r  which  is  equivalent  to  the  total 
transmittance,  thus 


0,,*)  Trl(z,9,)dr  .  (3.1)  st{z,)Tri(z,Bs )  -  ,«(«>) T^z,9s)  -  ,*(:)  .  (3.7) 


The  path  function  L.(z,9„0)  can  be  expressed  in 
terms  of  a  sun  component,  a  sky-terrain  component  and 
an  emitted  component 

L,(z,9s,d>)  -  }i(z)o-(z,0 ) 

+  J  L(z,6',4'Mz,/3')dn+Ltg(z)  ,  (3.2) 


LtU,9},<t>)  -  f jt) (tr(z,^l) 

+  j  Hz,9\<t>')ir{z,P')dn/s*{:) 

+  L.a (*)/,«<*))  .  (3.3) 

Let  the  last  two  terms  be  designated  by  the  function  C(z) 
such  that 

C(z)  -  A(*,*,.dV(a,0,)dft/J«(*)  +  £,a (*)/,«(*)  .  (3.4) 

Pubstituting  Eqs.  (3.3)  and  (3.4)  into  Eq.  (3.1)  we  get 
Ll(zSv<)>)  -  f  Jt(z,)T/.l(z,9J)[<r(z/,f3)  +  C(z,)]dr  .  (3.5) 


Now  substituting  into  Eq.  (3.S),  the  sun  irradiance  can  be 
taken  out  of  the  integral  and  we  get 

L^(z,9rd>)/Si(z)  -  f  1(<t(z,0)  +  C(z)]dr .  (3.8) 

M 

The  quantity  in  the  square  brackets  is  a  function  of 
altitude  and  not  of  zenith  angle,  therefore 

J*  [<r(z./3)  +  C(z)]dr 
00 

-  (m*(r,0,)/m„(z,O)]  J"|tr(z,j3)  +  C(j))<fe  .  (3.9) 

00 

where  ms.(z,0J)/m»(z,O)  is  the  relative  optica!  airmass  at 
9,.  For  altitudes  up  to  6  kilometers  we  can  use  the  sea 
level  value  for  sun  zenith  angles  0  to  86  degrees  so  for 
convenience  we  will  shorten  this  to  m(9,).  Substituting 
back  into  Eq.  (3.8)  we  get 

J#(z)m(0J)l  -  [cr (z ,/5)  +  C(z)]uz  .  (3.10) 

The  right  hand  term  can  be  separated  into  two  parts 

f  lo-iz.p)  +  C(z)]dz  -  J  <r(z,fi)dz  +  t  C(z)dz  .  (3.11) 


Now  let  us  define  an  optical  thickness  function  t„(z,0) 

such  that 


T„(i,p)  -  J"  rr <z,/£ 


The  integral  of  the  optical  thickness  function  over  An  is 
the  optical  thickness  due  to  scattering  ,/«(z) 


tM(r)  -  J* 


rm(it0)d(l 


checking  the  optical  stability.  Equation  (3. IS)  can  also  be 
written  as 


-  TJ:.0)m{ I.(~)|T00i2.fl)  +  /  C(z)dr)  .  (3.18) 

M 

Taking  the  log  of  both  sides  we  get 
loglZ,i(z,eJ,3)/m(eJ)}  -  m(»s)  logr»(z.O)  +  log  A  (3.19) 


since  the  optical  thickness  due  to  scattering  is  the  integral 
of  the  total  scattering  coefficient  with  altitude 


Substituting  Eq.  (3.12)  into  Eq.  (3.10)  we  get 


L*m (*,«,.*)/(,«(*)«(«,)] -rJtS)  +  f  C(:)dz  .  (3.1 


This  is  the  basic  equation  for  the  radiance  in  the 
solar  almucantar  and  how  it  relates  to  the  sun  scalar  irradi- 
ance,  the  optical  airiness,  the  optical  thickness  function 
and  one  additive  component  defined  by  Eq.  (3.4). 

3.2  Optical  Stability 


A  -,•(«>) Ir,,{z,/3)  +  J*  C(i)dz]  .  (3.! 


For  an  optically  stable  day,  A  should  be  constant  for 
a  given  angle  from  sun.  Thus  the  sky  radiance  at  a  con¬ 
stant  0  in  the  solar  almucantar  can  be  used  to  obtain  the 
beam  transmittance  when  measured  over  a  large  range  of 
airmass  values  on  an  optically  stable  day.  A  semi-log 
graph  of  ioglLl  [z,9,,0)/m{9,)]  versus  airmass  would  be 
linear  with  a  slope  of  log  transmittance  and  an  intercept  at 
log/4 .  The  transmittance  obtained  from  the  sun  radiance 
using  Eq.  (3.17)  and  the  transmittance  obtained  from  the 
sky  radiance  in  the  solar  almucantar  are  then  averaged  to 
obtain  the  most  accurate  value. 

Tashenov  (1970)  used  this  method  with  measure¬ 
ments  of  the  solar  aureole  to  obtain  spectral  transmittance 
in  the  region  410  to  733  nanometers. 


The  atmosphere  is  considered  to  be  optically  stable 
when  the  atmospheric  transmittance  does  not  change  with 
time.  When  the  day  is  optically  stable,  measurements  of 
the  apparent  sun  radiance  ,Lm(z,9„ 0)  made  over  a  large 
range  of  sun  zenith  angle  can  yield  a  good  measure  of 
total  transmittance  Tm(z, 0)  and  the  inherent  sun  radiance 
iLoM 

,Lm(z,9,.0)  -  ,LeMTm(t, 0)  .  (3.16) 

Taking  the  log  of  both  sides  of  Eq.  (3.16)  we  get 


Tests  to  Determine  Optical  Stability 

Pyaskovskaya-Fesenkova  (1970)  outlines  two  tests 
for  optical  stability  as  well  as  a  unique  equation  for 
transmittance  which  is  valid  only  for  optically  stable  days. 

The  first  test  is  to  graph  for  a  constant  scattering 
angle  the  ratio  I^(z,0,,/))/l,«(z)m(0,))  versus  m(0,). 
As  can  be  seen  from  Eq.  (3.15),  if  the  day  is  optically 
stable  the  right  hand  term  is  constant.  Therefore  this 
graph  should  result  in  a  horizontal  straight  line.  That  is, 
for  an  optically  stable  day 

Z.2(z,0r#)/[J«(z)m(0,)]  "  B  (3.21) 


g 


log  ,£„(*, 9,,0)  ■■  log  j Lq  (°°)  +  m(0,)  log  Tmiz, 0)  .  (3.17) 

Graphs  of  apparent  sun  radiance  on  semi  log  paper  of  log 
apparent  sun  radiance  versus  airmass  yield  the  typical 
Langley  (or  Bouguer)  graph  where  log  inherent  sun  radi¬ 
ance  is  the  intercept  and  log  transmittance  is  the  slope. 

Variations  of  Eq,  (3.15)  give  both  another  method 
for  obtaining  the  total  transmittance  and  several  ways  of 


B  -  r„,(z,/8) 


C(z)d! 


The  second  test  for  optical  stability  is  to  graph  the 
ratio  L»(z,0„/3)/,«(z)  versus  relative  airmass  for  a  con¬ 
stant  0.  The  graph  should  result  in  a  straight  line  going 
through  the  origin  on  an  optically  stable  day.  To  under¬ 
stand  this.  Eq,  (3.15)  is  rewritten  In  ihe  form 


lL  (z.#5.0)/s.(z)  -  m(«s)  5  .  (3.23) 


The  B  is  now  the  slope  and  the  intercept  is  zero. 

If  the  day  is  optically  stable  a  graph  of  solar  almu- 
cantar  radiance  at  a  constant  scattering  angle  versus  air- 
mass  will  yield  a  value  of  transmittance  as  Follows.  If  we 
write  Eq.  (3.21)  In  the  form 

L*m  {i,9,,fi)  -  ,«M  Tm(i. 0)m(s,)  m{»,)  B  (3.24) 

and  differentiate  (only  the  sky  radiance  and  the  air  mass 
are  variables,  the  inherent  sun  scalar  irradiance,  B  and  the 
transmittance  are  constant),  we  get 

dL*m{i,9,4 8) 

-  B  ,«(«*)  r„(*. or***’  d  m(0j)tl+m(01)/»7'.,(:,O))  .  (3.25) 

Equation  (3.2S)  is  only  valid  when  the  angle  from  sun  Is 
constant,  the  vertical  transmittance  from  space  to  sensor  Is 
constant  and  the  day  is  optically  stable.  Setting 
dLi  (z,0i,jS)  equal  to  zero  and  rearranging  we  get 

In T. (s,0)  *•  -l/m(0j)  (3.26) 

where  is  the  sun  angle  where  Z.i(z,0„/9)  is  a  max¬ 
imum.  Therefore  a  graph  of  solar  almucantar  radiance  at 
a  constant  angle  from  the  sun  versus  alrmass  will  Increase 
and  then  decrease  with  the  maximum  at  m(0j)  when  the 
day  is  optically  stable.  If  tho  curve  is  reasonably  well 
defined  with  sufficient  alrmass  values,  a  reasonable 
transmittance  value  can  be  obtained.  It  can  also  be  used 
as  a  third  teat  for  optical  stability.  The  values  of  transmit- 
tanoe  obtained  by  means  of  Eqs.  (3,17),  (3.19)  and  (3.26) 
ahould  be  in  agreement  on  an  optically  stable  day. 

Pyaskovskaya-Fesenkova  (1970)  describes  an  instru¬ 
ment  for  measuring  the  quantities  in  the  preceding  equa¬ 
tions.  It  la  essentially  a  radiance  photometer  with  a  three 
degree  held  of  view  for  the  solar  measurement.  For  the 
aureole  measurement  a  two  degree  to  three  degree  circular 
zone  Is  used.  Hor  article  also  includes  a  theoretical  test  of 
the  error  in  obtaining  transmittance  from  the  sun  radiance 
when  the  aureole  in  the  three  degree  Held  is  necessarily  a 
part  of  the  sun  radiance  measurement.  She  concluded  the 
error  to  be  negligible. 

3.3  Method  of  Determining  Spurious  Sun  Reflectance 

Pyaskovskaya-Fesenkova  (1970)  also  indicates  that 
the  two  methods  of  checking  for  optical  stability  will 
falsely  Indicate  Instability  If  reflections  in  the  radiometer 
telescope  are  not  negligible. 

Conversely,  once  It  is  established  that  the  day  is  opt¬ 
ically  stable,  the  above  tests  for  optical  stability  can  be 
used  to  test  whether  a  sky  radiance  measurement  contains 


spurious  sun  reflection.  Sky  radiances  near  the  sun  are 
difficult  to  measure  unless  the  sun  can  be  occulted.  For 
example,  recent  tests  made  with  the  C-130  sky  scanning 
photometer  indicate  spurious  sun  reflections  in  the  meas¬ 
urements  at  scattering  angles  0  to  25  degrees  from  the 
sun. 

To  illustrate  the  use  of  the  above  tests  for  optical 
stability  and  to  test  for  spurious  sun  reflections  in  a  sky 
radiance  photometer,  we  will  use  a  set  of  data  containing 
measurements  with  a  solar  transmissometer,  a  sky  scanner 
and  an  irradiometer.  These  measurements  were  made 
from  January  to  September  1964  from  a  platform  on  the 
roof  of  one  of  the  Visibility  Laboratory  buildings  at  Point 
Loma.  The  solar  transmissometer  measured  both  the 
center  3un  radiance  and  an  aureole  radiance  0.573  degrees 
from  the  sun  center.  The  general  procedure  was  to  make 
a  set  of  measurements  at  each  10  degree  increment  of  sun 
zenith  angle  and  at  noon. 

The  photometers  were  fitted  with  four  optical  filters 
so  that  there  were  two  narrow  flltered-sensor  passbands 
and  two  broad  passbands.  The  two  narrow  passbands  had 
mean  wavelengths  of  459  and  661  nanometers,  the  photo- 
pic  passband  had  a  mean  wavelength  of  560  nanometers, 
and  a  fourth  broad  passband  represented  an  unflltered 
sensor  with  a  mean  wavelength  of  505  nanometers.  Rela¬ 
tive  spectral  response  curves  for  these  filters  are  given  in 
Appendix  A. 

A  graph  of  the  apparent  center  sun  radiance  as  a 
function  of  relative  airmass  for  the  afternoon  of  2  Sep¬ 
tember  1964  is  given  in  Fig.  3-1.  The  straight  lines  for 
each  Alter  are  the  result  of  least  squares  flts  to  the  data. 
The  correlation  coefficients  were  high,  equal  to  or  greater 
than  0.99  and  the  day  Is  apparently  optically  stable.  The 
least  squares  transmittances  are  noted  on  the  graph. 


Fig.  3-1.  Langley  graph  of  measured  apparent  center  sun  radiance 
for  2  September  1964. 
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Solir  Almucantar  at  55  degrees 

In  order  to  be  certain  to  have  an  almucantar  sky 
radiance  measurement  without  sun  reflection  problems,  an 
angle  from  sun  of  5S  degrees  was  chosen.  The  sky 
scanner  data  are  for  a  fixed  set  of  zenith  angles  and 
azimuths  (A0  and  A0  “  S.62S  degrees).  Therefore  the 
almucantar  radiance  used  is  from  interpolations  to  the 
zenith  angle  of  the  sun  and  to  the  appropriate  azimuth  for 
ft  equal  to  55  degrees.  The  example  Is  done  for  only 
Filter  1  mean  wavelength  459  nanometers.  The 
almucantar/ m(9,)  graph  versus  airmsss  (Eq.  (3.19))  is 
given  in  Fig.  3-2,  The  transmittance  from  the  almucantar 
is  0.7S6  compared  to  the  transmittance  from  the  center 
sun  radiance  of  0.738  thus  there  is  only  a  2.4  percent 
difference.  The  average  0.747  is  considered  the  best 
measure  of  the  transmittance  if  the  day  is  optically  stable. 


n«.  3*2.  Almucantar  radiance  to  airmais  ratio  vs  relative  optical 
airman  for  2  September  1964. 


Now  for  the  three  tests  for  optical  stability. 
Although  we  do  not  have  a  measure  of  the  solar  scalar 
irradiance  ,«(0),  we  do  have  a  center  sun  radiance 
i £«.(0,9i,0)  which  is  related  to  the  scalar  irradiance  by 


f«(0)  -  1i..(o,el,o)Z)(3t)n1 .  (3.27) 

where  D(k)  is  the  center  to  average  sun  radiance  conver¬ 
sion  or  limb  darkening  factor  for  the  Alter  (see 
Appendix  A)  and  n,  it  the  solid  angle  subtended  by  the 
tun  for  2  September  1964.  If  we  substitute  the  center  sun 
radiance  for  the  solar  scalar  Irradiance  in  Eq.  (3.21)  we  get 


£*  (x,0l,/9)/{|Z.N(O,0l,O)m(0J)]  -  B' ,  (3.28) 


where 


Since  D(\)  is  a  constant  for  a  filter  and  n,  is  constant  for 
any  given  day,  the  two  tests  can  now  be  made  using 
Eqs.  (3.28)  and  (3.30)  instead  of  Eqs.  (3.21)  and  (3.23). 

Fig.  3-3  contains  the  appropriate  graphs  for  the  two 
tests.  Although  not  perfect,  the  data  in  Fig.  3- 3a  is  rea¬ 
sonably  represented  by  a  horizontal  straight  line.  Simi¬ 
larly,  the  data  in  Fig.  3-3b  are  reasonably  represented  by  a 
straight  line  going  through  the  origin. 

oa. 
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Fig.  3-3a.  Oraphi  to  teat  atmospheric  optical  stability  for  2  Sep¬ 
tember  1964,  filter  1  (A  -  459  nm),  Test  One,  Eq.(3.28). 


- 1 - 1 - r- 

- 1 - 1 - 

TEST  ONE  . 

K  -  459nm 

*  EQ.  (3.28) 

- 

• 

_ i  ■  ‘ 

_ I _ 1 

RELATIVE  OPTICAL  AIRMASS,  m  ( 0, ) 

Fig.  3-3b.  Graphs  to  test  atmospheric  optical  stability  Tor  2  Sep¬ 
tember  1964,  Alter  1  (A  -  459  am),  Test  Two,  Eq.(3.30). 


B'-D(K)n,B.  (3.29) 

Similarly  Eq.  (3.23)  becomes 

Li  (x,0J,/9)/J/.„(r,0J,O)  -  miS,)  B'  .  (3.30) 


The  third  test  is  a  plot  of  the  almucantar  sky  radi¬ 
ance  graphed  against  expl-l/m(0,)]  from  Eq.  (3.26)  as 
shown  in  Fig.  3-4.  The  maximum  sky  radiance  should  lie 
at  a  transmittance  of  approximately  0.75.  Although  the 
noon  value  of  sky  radiance  is  inconsistent  with  the  rest 
(hooks  upward  near  air  mass  one),  and  there  are  no  mcas- 
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ured  iky  radiances  at  the  exact  the  sky  radiances  can 
be  reasonably  represented  by  a  curve  which  maximizes  at 
a  transmittance  of  0.75  (dashed  line). 

120 1 - 1 - . - . - 1 - 1 


At  i  I...  . I  - . I  —.—...I 

u0  .20  .40  ,60  80  1.00 

VERTICAL  TRANSMITTANCE 
T.(0, 0)  -  EXP|-l/m(e,)l 

Fla.  3*4*  Almucantar  sky  radiance  at  55*  vs  vertical  transmittance 
(Eq.  3.26)  for  2  September  1964,  filter  1  (k  -  459  nm). 

Thus  we  conclude  that  the  atmosphere  on  2  Sep* 
tember  1964  In  the  vicinity  of  the  Visibility  Laboratory 
was  reasonably  stable  optically.  Also  we  conclude  that 
scanner  radiances  may  be  used  to  obtain  almucantar 
values  by  Interpolation  for  use  in  the  preceding  equations. 

Now  that  it  is  established  that  the  day  is  optically 
stable,  wo  will  use  these  same  equations  to  test  for  sun 
reflections  in  the  aureole  photometer. 

AuartiBiHinH 

The  solar  aureole  radiance  was  measured  0.573 
detrees  from  the  center  of  the  sun  without  occulting  the 
sun.  A  graph  of  the  ratio  of  the  apparent  aureole  radiance 
to  the  relative  airmass  versus  airmass  Is  given  in  Fig.  3-5. 


fig,  3*6a.  Graph*  to  test  for  spurious  sun  reflections  in  eureole 
radiance  measurements,  2  September  1964,  Test  One,  Eq.(3.28). 


Tig.  3*5.  Aureole  radiance  to  airmass  ratio  vi  relative  optica,  air¬ 
mass  for  2  September  1964. 


The  data  for  the  three  Niters  look  reasonably  linear  and 
the  correlation  coefficients  are  0.99.  However,  the  least 
squares  transmittances,  which  are  noted  on  the  graph, 
vary  markedly  from  the  values  derived  from  the  sun  radi¬ 
ance  and  the  almucantar  at  55  degrees  scatterinlg  angle. 

The  two  tests  for  spurious  sun  reflection  using 
Eqi.  (3.28)  and  (3.30),  fall  miserably,  see  Fig.  3-6a  and  b. 
In  Fig.  3*6a  the  B'  increases  dramatically  as  airmass  gets 
near  one  instead  of  staying  constant.  In  Fig.  3-6b  the  data 
can  be  fit  to  straight  lines  but  they  clearly  do  not  go 
through  the  origin. 


Fig,  3-6b.  Graph*  to  lei,  for  ipurlous  tun  reflections  in  eureole 
radiance  measurements,  2  September  1964,  Ttet  Two,  Eq.(3,30). 
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The  final  teat  it  the  graph  of  the  solar  aureole  radi¬ 
ance  as  a  function  of  expt-l/mW,)]  as  shown  in  Fig.  3-7. 
The  aureole  peaks  at  transmittances  still  lower  than  the 
values  derived  using  the  aureole  and  noted  in  Fig.  3-5. 
The  transmittances  from  Eqs.  (3.17),  (3.19)  and  (3.26) 
are  clearly  not  in  agreement.  Thus  we  conclude  that  the 
solar  aureole  measurements  include  a  spurious  sun 
reflection  component  In  all  three  Alters. 


The  last  term  in  Eq.  (3.32)  can  be  rewritten  substi¬ 
tuting  in  Eq.  (3.4)  for  COr) 

ifcMMn  -III  lLU.0\ii')(T{t,(i')/st(2)]dttdz  dl\ 
+  jf/  (£% <!)/,«(*)]*  dn  .  (3.33) 


IF  TM(Q,0)  -  EXP  I— 1/ni  (•,)  1 

fig.  VJ.  Solar  aureole  radiance  vs  vertical  transmittance  (Eq. 
3.26)  for  2  September  1964. 

Any  iky  radiance  measurement  made  near  the  sun 
without  occultation  of  the  sun  Is  questionable  and  should 
be  tested  in  some  fashion.  The  above  method  appears  to 
be  a  valid  testing  procedure. 

This  test  can  be  run  for  the  sky  scanner  used  on  2 
September  1964  for  various  angles  from  the  sun  smaller 
than  S3  degrees  to  establish  which  scattering  angles  have 
valid  sky  radiance  measurements  and  which  contain  spuri¬ 
ous  sun  reflections. 


3.4  Large  Sen  Zenith  Angle 

When  the  sun  zenith  angle  is  relatively  large,  the 
solar  slmucantar  contains  a  large  enough  range  of  scatter¬ 
ing  angles  so  that  Eq.  (3.13)  can  be  Integrated  over  4n 
thus 


2*/  Lm  (*.•,. /J)sln/M0/l  1«(x)m(SJ)] 
0 

C(:)dzdCl  . 

Substituting  in  Eq.  (3.13)  we  get 

2#/  L fi)iin0d0/ 1 

9 

-,/„(!)  +  j  f  cu)didn  . 


(3.31) 


(3.32) 


This  integrates  to 

f  J  C(z)dt  dll  -  f  {U(f)-,((f)]i(f )/,•(*))  dz 

T#  ee  ee 

+  J  [*.(*)/.«<*»*  .  (3.34) 

m  a 

But  the  first  term  in  Eq.  (3.34)  can  also  be  written  as 


J  (l»(r)-J»(r)|j(r)/,«(r)}  dr - 
/ l«(i)r(r)/,*(r)]dr  -  ,lm(*)  ■  (3.35) 

«n 

Therefore  combining  Eqs.  (3.32),  (3.34)  and  (3.35)  we  get 

J  Z.I(;r,SJ,0)»ln0d/}/l1«(»)m(0J)) 

0 

-/  («(*)*<*)/,«(*))*+  f|».(z)/ .1(1)]*  .  (3.36) 

"  ee  0 

By  multiplying  both  sides  by  the  sun  scalar  irradiance  out 
of  the  atmosphere  ,<(«>),  this  can  also  be  written  as 


2ir  J"  Lm(r,6,,0)tin0dfi/lrm(i,«t)m(9t)] 

9 

■|  exp(/  «(r)rfr)li(r)j(r)  +  «,  (z)}dk  .  (3.37) 

ee  ee  ® 


No  Absorption 

When  there  is  no  absorption  the  total  scalar  Irradi¬ 
ance  «(z)  does  not  change  with  altitude,  therefore  it  may 
be  taken  out  of  the  integral.  Also  the  attenuation 
coefficient  a(z)  is  equal  to  the  scattering  coefficient  j(z). 
In  addition  there  is  no  emission  thus  m.d(r)  equals  zero. 
Making  these  substitutions  into  Eq.  (3.37)  and  transposing 
m(9,)  we  have 


2ir 


J 


-«(;) 


s(i)dr]s(:)m(tis)d:  . 


(3.38) 
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For  the  sun  angles  where  m(6,)  equals  stc9,  and  dr 
equals  aec9  dx ,  wp  now  have  an  integral  of  the  form  e“du 
where  u  equals  j  s(s)sec9,dz  and  du  equals  s(z)sec9,dz. 
Therefore,  we  need  a  large  sun  zenith  angle  to  get  full 
coverage  of  the  scattering  angles  0  to  180  decrees  but  we 
need  the  sun  zenith  angle  to  be  equal  to  or  less  than  70 
degrees  so  that  .the  airtMss  is  equal  to  sec9,.  These  two 
requirements  are  a  bit  contradictory.  Assuming  the  two 
requirements  reasonably  compatible  at  70  to  80  degrees 
senlth  angle,  we  can  integrate  to  get 

2w]‘  Lm  (»,*,, 0)*in0  dfi/T„(i,9,) 

9 

-S(*)lr.(», g,)-1-!].  (3.39) 


Transposing  the  transmittance  terms  and  simplifying  we 
have 


U I  C<*,al^)slnflrfg/(l-r.(i, •,))«•<*)  .  (3.39) 

Whan  there  is  no  abeorption  the  normalized  equili¬ 
brium  radiance  from  Eq.  (2.49)  is 

£  y*,e,d)rffl  -  a(s) .  (3.40) 


nierefon,  Eqa.  (3.39)  and  (3.40)  Imply  that  for  the  solar 
almucantar  at  relatively  large  sun  zenith  angles, 

V*,#,,*)  -  t:(«,#|1/i)/(i-r.(i,sl)l .  (3.41) 


4.  VOLUME  SCATTERING  FUNCTION 


4,1  Integrating  Naghalametar 

The  Visibility  Laboratory  integrating  nephelometer 
was  first  deecrlbed  in  Duntiey  era/.  (1970).  Measure¬ 
ments  were  made  on  the  ground  and  airborne  by  the  Visi¬ 
bility  Laboratory  with  this  instrument  from  1968  through 
1973  in  locations  in  Thailand,  the  continental  United 
Statea  and  in  Europe.  From  1977  through  1980  similar 
measurements  were  made  in  Europe  with  a  comparable 
folded  path  instrument  described  In  Duntiey  *t  ai  (1973). 

The  Integrating  nophelometer  measures  three 
scattering  properties:  the  total  volume  scattering  coefficient 
j(i),  and  the  volume  scattering  function  9)  at  30  and 
130  degress  scattering  angle.  Measurements  between 
1968  and  1910  were  made  with  various  narrow  and  broad 
band  speotral  sensitivities  in  the  visible  spectrum  and  very 
near  infrared,  with  mean  waveienths  from  478  to  763 
nanomoters. 


When  the  integrating  nephelometer  data  are  put  into 
the  form  of  the  ratio  to  the  Rayleigh  scattering,  some  of 
the  differences  due  to  wavelength  and  air  density  changes 
with  altitude  disappear.  Let  us  define  a  ratio  Q(z)  as  the 
ratio  of  the  total  volume  scattering  coefficient  to  the  Ray¬ 
leigh  volume  scattering  coefficient  *i(z). 

Q(z)  -  i(r)/*i(r)  .  (4.1) 

Similarly  let  us  define  a  ratio  3(z,/3)  as  the  ratio  of  the 
total  volume  scattering  function  to  the  Rayleigh  volume 
scattering  function  «cr(z,/9) 

2<i,0)  ■*  o-(i,0)/*tr(*,0)  .  (4.2) 

When  the  scattering  function  ratio  2(z,/9)  is  graphed 
as  a  function  of  the  total  scattering  coefficient  ratio  <?(*), 
a  single  function  results  for  each  scattering  angle  /9, 
regardless  of  the  Alter,  or  whether  the  data  are  ground 
level  as  in  Fig.  4-1  or  airborne  u  in  Fig.  4-2.  All  the 
Integrating  nephelometer  data  from  1968  through  1980  are 
similar  regardless  of  location  of  measurement  or  filter. 

The  superimposed  curves  are  from  the  catalog  of 
normalized  volume  scattering  functions  from  Bartenova 
(1960). 


Fig.  4-1.  Comparison  of  multl-ipectral  volume  scattering  function 
ratios  measured  by  Visibility  Laboratory  ground-based  nephelom- 
alar  with  the  photopic  ratioa  from  Barteneva  (I960). 
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OPTICAL  SCATTERING  RATIO,  Q(t) 

Pig.  44.  Companion  of  multi-spectral  volume  Kettering  function 
rtlloi  meuured  by  Viilbillty  Laboratory  airborne  nephelometer 
with  the  photoplc  ratloa  from  Barteneva  <1960). 


4.2  Barteneva  Catalog 

Barteneva  (1960)  provided  a  catalog  of  normalized 
directional  acattering  functions  <r(fi)h  for  the  photoplc 
aenaor  baaed  upon  624  measurements  made  with  a 
nephelometer  during  19SS  through  19S8  In  various  loca¬ 
tions  in  the  USSR  and  at  sea.  For  each  volume  scattering 
function  she  givee  a  range  of  total  volume  scattering 
coefficient  appropriate  to  that  function. 

The  range  of  total  volume  scattering  coefficients  Is 
from  Rayleigh  to  heavy  fog,  for  ground  level.  Her  catalog 
of  gradual  (Unctions,  classes  I  through  10,  are  appropriate 
for  very  clear  to  heavy  haze,  but  not  fog.  This  range  of 
functions  Is  applicable  to  the  troposphere  under  conditions 
appropriate  for  flying  using  visual  flight  rules.  This  also 
describes  the  conditions  under  which  the  aircraft  flew  dur¬ 
ing  the  Visibility  Laboratory  data  missions.  Since  the 
ground-based  data  were  taken  in  coi\junction  with  the 
flights,  those  same  conditions  are  also  appropriate  for  the 
ground-based  data. 

The  Barteneva  catalog  was  first  interpreted  in  terms 
of  the  median  total  volume  scattering  value  for  each 
scattering  function  clau.  Then  the  ratio  Q  was  computed 
for  these  median  a  values.  Similarly,  the  value  of  the 
median  volume  scattering  function  <r(0,/9)  for  30  and  ISO 
degrees  was  computed  by  multiplying  the  normalized  func¬ 
tions  by  the  median  a  values.  Then  the  ratio  to  Rayleigh 


V(z,/3)  was  computed  and  the  curves  in  Figs.  4-1  und  4-2 
graphed  for  the  Barteneva  data. 

All  the  Visibility  Laboratory  nephelometer  data  for 
total  volume  scattering  and  the  scattering  functions  at  30 
und  ISO  degrees  agree  with  the  Burtencvu  catalog  of  gra¬ 
dual  functions.  Thus  it  has  been  assumed  that  the  Bar¬ 
teneva  catalog  is  appropriate  for  specifying  the  directional 
scattering  for  all  the  scattering  angles  0  to  180  degrees  for 
various  sensors  in  the  visible  spectrum  in  the  troposphere. 

The  further  step  of  representing  the  catalog  by  equa¬ 
tions  has  been  made  by  Johnson  and  Hering  (1981a).  A 
full  analytic  representation  in  the  form  of  Henyey- 
Greenstein  functions  is  presented  in  Hering  (1981b). 

Graphs  of  the  normalized  volume  scattering  func¬ 
tions  as  a  function  of  scattering  angle  for  the  Barteneva 
catalog  are  given  in  Fig.  4-3  for  the  gradual  classes  and 
Fig.  4-4  for  the  steeper  classes. 


SCATTERING  ANGLE,  0 

Fig,  4-3.  Normalized  volume  scattering  function  for  gradual 
classes  horn  Fig.  I  Barteneva  (I960). 


30  90  ISO 

SCATTERING  ANOLE.  0 


Fig.  4-4.  Normalized  volume  scattering  function  far  steeper 
classes  from  Fig.  2  Barteneva  (I960). 


(4.7) 


4.3  Avtragt  Volume  Scattering  Function  at  55° 

The  sverage  normalized  volume  scattering  function 

is 


|<r(j3)/j)dn/  (  dn-  l/(4ir)  -  .08  .  (4.3) 


The  Raylalgh  normalized  volume  scattering  function  can 
be  computed  from  the  approximation 


ft<r(/S)/*i  -  3(1  +  coi,/3)/(l6«r)  .  (4.4) 

Equations  (4.3)  and  (4.4)  were  combined  and  the  ensuing 
aquation  solved  to  obtain  the  scattering  angle  at  which  the 
Rayleigh  normalized  volume  scattering  function  is  equal  to 
the  average  value  of  l/4ir.  It  waa  found  to  be  53  degrees 


*<r(5J )/#s  -  1/(4#)  .  (4.5) 


The  graphs  in  Fip.  4-3  and  4*4  are  marked  at 
scattering  angle  35  degrees  and  at  the  normalized  volume 
scattering  function  equal  to  0.08.  From  these  graphs  we 
see  that  the  scattering  angle  of  S3  degrees  yields  a  reason¬ 
able  approximation  of  the  average  volume  scattering  (Unc¬ 
tion  for  all  the  scattering  (Unctions  except  the  gradual 
daaeea  9  and  10,  that  is, 


<r(SS)/f  -  1/(4#)  .  (4.6) 

In  reviewing  the  graphs  of  measured  nephalometer 
data,  It  is  noted  that  none  of  the  total  scattering  coefficient 
ratios  (?(»)  exceeded  100.  That  is,  gradual  classes  1 
through  I  (Q(8)  ■  114)  were  encountered  but  never  9  or 
10  (Q(9)  "  201  and  Q(10) »  237),  Therefore  for  model¬ 
ing  the  atmoepheres  measured  from  1968  through  1981  in 
the  visible  spectrum,  It  it  reasonable  to  use  Eq.  (4.6)  as 
part  of  the  model  assumptions. 


4.4  Mat  bed  of  Obtaining  Single  Scattering  Albedo 
from  Horison  Sky 

Roessier  and  Paxvog  (1981)  developed  equations 
for  the  horizon  aky  luminance  as  a  function  of  the  single 
scattering  albedo  after  assuming  both  an  Isotropic  scatter¬ 
ing  (Unction  and  an  isotropic  luminance  distribution.  The 
technique  can  be  made  more  generally  applicable  by 
developing  the  equations  for  the  horizon  it  a  scattering 
angle  of  33  degress  from  the  sun. 


Hgflwn  Sky  n  Ml  torn  Sun 

The  equation  for  the  path  function  at  55  degrees 
scattering  angle  can  be  expressed  is  a  (Unction  of  the  sun 
component  separate  from  the  sky-terrain  component  in  a 
manner  similar  to  Eq.  (3.2) 


55  (1)  »  ,t (r)cr(r,S5) 

+  f  [2,(1') dll  +  L.g(2)  . 


Substituting  Eq.  (4.6)  into  Eq.  (4.7)  and  assuming  the 
sky-terrain  or  diffuse  scalar  irradiance  rfi(z)  Is  Isotropic  we 
get 

£.«(?, 0, 55  fi)  «■  {c(t)t(z)/(4fr) 

+  rf«(z)i(z)/(4#>  +  £.,fl(z).  (4.8) 

Since  the  total  scalar  irradiance  is  the  sun  plus  the  sky- 
terrain  scalar  Irradiance,  Eq.  (4.8)  can  be  written  as 

£.•(1,9.35  fi)  -  •(*)*<!)/( 4#)  +  L%  <*)  .  (4.9) 

0 

Dividing  both  sides  by  the  attenuation  coefficient  and 
expressing  the  emitted  path  function  according  to 
Eq.  (2.7)  we  get  an  equation  for  the  equilibrium  radiance 

£.,  (*,0,55  0)  ■  «(*)t(*)/(a(r)4#] 

+  s(i)lU,D/e(i).  (4.10) 


This  can  be  expressed  In  terms  of  the  single  scattering 
albedo  a*  (from  Eq.  (2.39))  as 


£.,(1,9,55  fi)  -  w(i)«(i)/(4ir) 

+  li-w(*)U<\,r) ,  (4.U) 

The  general  equation  for  the  horizon  sky  radiance  at 
55  degrees  scattering  angle  from  Eq.  (2.57)  is 

£« (#.90,55)3)  -  £.,(*, 90,55/8)[l~r..(*.90)l  .  (4.12) 

Therefore  from  Eq.  (4.11)  the  horizon  sky  radiance  at  55 

degrees  is 

il  (1,90,55/8) 

-  |w(*)*(*)/(4#)  +  (l-«U))£.<X.r))Il-r.(xl90)l  .  (4.13) 

The  above  equation  assumes  that  the  attenuation 
coefficient  Is  horizontally  isotropic  and  therefore  the  hor¬ 
izon  sky  must  be  cloud  free. 


Visible  Spectrum 

In  the  visible  spectrum  LU.D  is  negligible.  Also 
for  the  photopic  at  sea  level  f»(0,90)  is  negligible  and 
Eq.  (4.13)  becomes 
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1*  (0.90,550)  -<«(0)«(0)/4w  . 


(4.14) 

Thus  the  single  scattering  albedo  w  can  be  determined  by 
measuring  the  horizon  sky  at  55  degrees  and  the  total 
scalar  irradiance  at  sea  level 

«(0)-i;<0,90,J3/8)4ir/«(0)  .  (4.15) 

Since  this  Is  an  approximative  method,  the  precision  level 
of  the  resultant  *»(0)  is  not  expected  to  be  high.  How¬ 
ever,  it  can  still  be  quite  useful  as  one  measure  of  the  sin¬ 
gle  scattering  albedo. 

When  there  Is  no  absorption  u  -  1  and  the  horizon 
sky  should  equal  the  average  scalar  irradiance  «(0)/(4ir) 


(0,90,550)  «■  «(0)/(4ir) .  (4.16) 


4,5  Method  of  Obtaining  Scattering  Transmittance 
hem  Sky  Radiance  Ration  at  55  Degrees 

When  abeorptlon  Is  negligible,  the  general  equation 
for  the  equilibrium  radlanoe  at  S3  degrees  scattering  angle 
Bq.  (4.11)  reduces  to 

1,  Or, 0,55  0)  -  «(i)/4ir  .  (4.17) 

Also  when  absorption  is  negligible  the  total  scalar  Irradl* 
ance  «(r)  is  constant  with  altitude,  therefore  the  equili¬ 
brium  radiance  at  S3  degrees  Is  also  constant  with  altitude. 
Rearranging  the  equilibrium  radiance  form  of  the  equation 
of  transfer  Eq.  (2.17)  we  get 

(£(r,0,55  0)-4,<i ,0,55  0)]~>  4L  (i ,0,55  0) 

--!(*)*,  (4.18) 

where  the  total  scattering  coefficient  s(x)  has  been  substi¬ 
tuted  for  the  attenuation  coefficient  since  absorption  Is 
negligible.  Since  the  equilibrium  radiance  is  a  constant 
Bq.  (4.11)  an  be  Integrated  with  respect  to  altitude 


|  (I (z.0.55  0)-i,  (*,0,53  0))"'  atlU.0,55  0) 

•* 

--/»(*)<*.  (4.19) 


The  result  of  the  Integration  is 

l,(t ,0,55  0)  -  £o(f,,0,55  0) ,  r,(*,0) 

+  1,(1,0,55  0)11  -//>,  0)1  (4.20) 


where ,  T,(z,0)  Is  the  transmittance  due  to  scattering. 


Sky  Radiance  Ratio  at  55° 

For  sky  radiance  at  55  degrees  scattering  angle,  the 
inherent  radiance  1„  is  equal  to  zero  and  Fq.  (4.20) 
becomes 

L*m (z.0,55  0)  -  !,(*, 4,55 0)|l  -jf.tr.fl)]  .  (4.21) 

Now  when  55  degree  sky  radiances  at  two  zenith 
angles  0  and  9'  are  ratioed,  since  the  equilibrium  radiance 
is  equivalent,  the  sky  radiance  ratio  is  solely  a  function  of 
the  scattering  transmittance 

1m  (z.0.55  0)/1m  (z.0',55  0) 

-  |l  -  ,  rjx.OFW)  /  [1  -  ,  rM^O)'*"®'1)  .  (4.22) 


Although  Eq.  (4.22)  cannot  be  solved  for  directly,  it  can 
be  solved  by  iterative  means.  Error  analysis  indicates  that 
the  zenith  angle  difference  0-0'  should  be  large  to 
minimize  the  error  in  the  resultant  vertical  transmittonce 
,r„(x,0). 

Visible  Spectrum 

In  the  visible  spectrum,  absorption  Is  negligible 
except  for  ozone.  The  total  transmittance  In  the  visible 
spectrum  would  thus  be  approximated  by  the  product  of 
the  satterlng  transmittance  times  the  ozone  transmittance 


r.(z, o)  - ,  r,(».o)  fl  r«(i,o) ,  (4.23) 

To  Illustrate  the  use  of  Eqs.  (4.22)  and  (4.23),  we 
again  refer  to  the  sky  radiance  and  solar  transmissometer 
data  taken  on  the  Visibility  Laboratory  rooftop  in  1964. 
The  atmosphere  on  2  September  1964  was  optically  stable 
during  the  afternoon  as  Illustrated  In  Section  3.3.  Thus 
evaluation  of  the  satterlng  transmittances  throughout  the 
afternoon  should  give  us  an  estimate  of  the  precision  of 
the  estimation  method  for  broad  band  sensors  in  the  visi¬ 
ble.  Equation  (4.22)  was  evaluated  for  all  sky  radiances  at 
55  degrees  plus  or  minus  2.5  degrees  satterlng  angle  for  0 
from  81.6  to  64.7  degrees  and  0‘  from  64.7  to  2.8  degrees. 
Later  error  analyses  indlated  that  some  of  these  zenith 
angle  combinations  are  less  error  prone  than  others. 
Hence,  the  averages  In  Table  4. 1  are  less  accurate  than  can 
be  obtained  with  a  smaller  yet  better  selection  of  zenith 
angle  combinations.  However,  they  are  presented  herein 
as  a  first  approximation. 

The  transmittance  ratios  in  Table  4.1  are  the 
transmittances  based  on  Eq.  (4.22)  times  the  ozone 
transmittance  divided  by  the  transmittance  from  the  solar 
transmissometer.  The  ozone  transmittance  is  noted  at  the 
bottom  of  each  column.  The  transmittance  comparison  is 
best  for  Filter  1  which  is  a  relatively  narrow  band  filter 
with  the  least  absorption.  The  near  noon  data  compare 
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Tablt  4.1.  Comparison  of  transmittance  determination  procedures 
(sky  radiance  ratio  method  vs  solar  transmissometer). 


Nomlnil  Sun 
Zenith  Anile 

TOTAL  transmittance  ratio 

(From  Eq  (4.22,  and  4  23)  grid  Solar  Trantmiuometcr  Measurements) 

(dHfMft) 

Filler  1 

i-  459** 

Filter  2 

A  -  505  nm 

Filler  ) 

A  -  560  nm 

Filter  4 

A  -  661  mi 

14 

.954 

.HI 

.194 

.900 

SO 

.944 

193 

.194 

.905 

40 

1.045 

1.011 

1.050 

1.0)1 

M 

1.045 

t.0!) 

1.041 

1.047 

to 

1.0)4 

1.004 

1.0)1 

1.044 

TO 

I.0M 

1.014 

1.04) 

1.050 

so 

100) 

.99(1 

1.0)1 

1.045 

Avert* 

toot 

,9TS 

.995 

1.004 

tumlvd 

Deviation 

045 

.057 

.099 

.070 

Omm 
ContflMIcn 
,Tm  (0.0) 

ttt 

.914 

.969 

.971 

tout  w« II.  The  eccuncy  ia  probably  auflklent  to  warrant 
further  development  of  the  method  to  improve  the  prod* 
aion  level. 

Uee  of  aky  radiancea  at  S3  degrees  scattering  angle 
should  be  an  Improvement  on  the  sky  ratio  method 
described  In  Duntley  etal.  (1972)  Section  2.1  and  In 
Duntley  ttal.  (1971)  Section  2.2.  That  method,  which 
stemmed  from  the  nomographic  method  of  Kushpil'  and 
Petrova  (1971),  used  aky  radianou  at  all  scattering  angles. 
Kushpil’  and  Petrova  (1971)  suggested  use  of  ratios  at 
S7.2  degrees  scattering  angle  for  the  visible  spectrum  and 
S3.9  degrees  for  the  near  infrared  portion  of  the  spectrum 
but  did  not  give  a  theoretical  buls  for  the  nomograph  or 
for  the  selection  of  these  angles. 


4.C  Met  bed  ef  Obtaining  Aereeol  Optical  Thickness 
hen  Solar  Alwucantar 

The  bulc  equation  for  the  sky  radiance  in  the  solar 
almucantar  Eq.  (3.15)  can  be  rewritten  expressing  the  opt¬ 
ical  thlckneu  function  in  two  components,  the  Rayleigh  or 
molecular  *r,(*,/ S)  and  the  Mlt  vr.(r,/3) 


S)  /  (,•<»)»»(*,)] 

«  *t„(i,/3)  +  ^t.(i,/})  +  J*  C(i)dr  .  (4.24) 

■a 


For  simplicity  the  left  hand  ratio  is  defined  as  a  function 


-  Li  0-,9,,/S)  /  .  (4.25) 

if  we  have  in  the  almucantar  a  large  enough  range 
of  scattering  angles,  both  sides  of  Eq.  (4.24)  can  be 
Integrated  over  4tr  and  Eq.  (3.13)  substituted  In  to  obtain 


J  m(s./3)  »lni8  dfi  -  arji)+  wr„(*) 

0 

+  2ff  J  J  C(r)dr  t\t\fidp  .  (4.26) 


where  *r»(jt)  is  the  Rayleigh  optical  thickness  and  uu(i) 
is  the  Mle  optical  thickness.  The  last  term  in  both 
Eqs.  (4.24)  and  (4.26)  is  the  diffuse  component  due  to  the 
contribution  of  the  relatively  non-directional  sky,  apparent 
terrain  radiance  and  emission.  All  the  equations  up  to  this 
point  are  not  approximations  and  hold  for  an  atmosphere 
with  or  without  absorption  and/or  emission.  Thus,  they 
are  equally  as  valid  in  the  infrared. 

Livshits  and  Pavlov  (1970)  make  the  assumption 
that  in  the  photopic  or  visible  portion  of  the  spectrum,  the 
sky  luminance  or  radiance  in  the  almucantar  for  scattering 
angles  90  to  180  degrees  is  essentially  a  function  of  the 
Rayleigh  optical  thickness  function  and  the  diffuse  com¬ 
ponent,  and  that  the  diffuse  component  is  the  same  for  0 
to  90  degrees  as  it  is  for  90  to  180  degrees  scattering 
angle.  This  assumes  homogeneity  In  the  almucantar  and 
thus  is  limited  to  cloudless  skies.  Thus  since  the  Rayleigh 
scattering  is  symmetrical  about  90  degrees 
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tLitiieUal 


2[2w  J  4(2,0)  sln/8</0] 
in 

-  +  2ir  J*  j"  Ciz)dz s!n0d/3  .  (4.27) 

0  o* 

Alio  it  follows  (hit  the  Mie  optical  thickness  is  the  Integral 
of  the  forward  portion  of  the  function  less  the  Integral  of 
the  back  portion 


M »„(*)  «•  2ir  J  n(isP)  sin/3  4/9 
0 

-2ir  J  m(i, /3)  sln/9 4/9  .  (4.28) 

in 


Livshits  and  Pavlov  (1970)  also  obtained  the  total 
optical  thickness  tm(x)  from 


r„(i)  -  a  »»(*)  +  huM  +  fl  tmh>  •  (4.29) 


whore  ,u(x)  is  tho  optical  thickness  due  to  absorption. 
They  obtain  the  Rayleigh  thickness  theoretically,  and 
assume  the  absorption  is  primarily  due  to  ozone  in  the 
visible  wavelengths.  Ozone  absorption  was  assumed  to  be 
the  yearly  average  value  for  a  given  wavelength.  They  got 
good  comparison  to  totsl  optical  thickness  values  obtained 
with  the  Langley  or  Bouger  method  for  short  wavelengths 
where  there  are  no  other  absorption  bands.  They  did  not 
compare  as  well  at  the  longer  wavelengths.  We  suspect 
absorption  other  than  for  ozone  In  the  longer  wavelengths 
In  the  visible  rather  than  an  error  in  the  aerosol  optical 
thickness. 


agmmnt.AMH 

Tho  disadvantage  of  Eq.  (4.28)  Is  the  need  for  a 
large  sun  zenith  angle  so  that  a  large  range  of  scattering 
angles  are  available  in  the  solar  almucantar  radiances. 

Since  the  ratio  of  the  volume  scattering  function  at 
33  degrees  to  the  scattering  coefficient  is  l/(4ir)  for  the 
Rayleigh  and  the  total  scattering  (Eqs.  (4.5)  and  (4.6)),  it 
must  also  be  l/(4tr)  for  Mie  scattering  as  well 


-  JI»(s1J5)/lfj(«)  -  l/(4tr)  .  (4.30) 

Thus  the  ratio  of  the  optical  thickness  function  at  53 
degrees  to  the  optical  thickness  must  also  be  l/(4ir)  for 
each  component  as  well  as  tite  total  scattering 

*».(*, J 5)/ttmh)  -  Tm(z,SS)/,fm(t) 

-  wr,(t,J5)/wfM(*)  -  l/(4tr)  .  (4.31) 


Now  we  can  express  the  Mie  optical  thickness  function  in 
terms  of  the  Mie  optical  thickness  by  rearranging 
Eq.  (4.31) 

wr„(2,55)  -  M/„(z)/(4ir)  .  (4.32) 

Using  the  same  rationale  as  Livshits  and  Pavlov 
(1970),  that  the  angles  90  to  180  degrees  in  scattering 
angle  describe  the  diffiise  and  Rayleigh  components  for 
the  angles  0  to  90  degrees,  we  can  say  that 


/u(x,180  -  55)  -  ftT„(r,55)  +  j*  C(t)dz  .  (4.33) 


Also  from  Eq.  (4.24)  we  have 

m(z,5S)  -  MTm(t,5S)  +  KTm(i,$S)  +  j"  C(i)dz  .  (4.34) 


Subtracting  Eq.  (4.33)  from  Eq.  (4.34)  we  get 

tfT-U.JS)  -  m(,,55)  -  4(1,125) .  (4.35) 

Now  combining  with  Eq.  (4.32)  we  get 

*,r..(*)  -  4»rU(r, 55)  -»*(«,  125))  ,  (4.36) 

The  solar  almucantar  contains  both  55  and  125 
degrees  scattering  angle  for  sun  zenith  angles  62.5  to  90 
degrees.  Thus  for  these  solar  zenith  angles  a  measure* 
ment  of  the  scalar  sun  irradiance  and  sky  radiances  in  thr 
solar  almucantar  at  55  and  125  degrees  from  the  sun 
would  yield  a  measure  of  the  Mie  optical  thickness 
through  the  atmosphere.  As  long  as  the  sky  radiance  dis¬ 
tribution  shows  brighter  areas  0  to  90  degrees  than  90  to 
180  degrees,  this  scheme  will  work.  Equation  (4.36)  also 
Is  valid  when  atmospheric  emission  is  present  as  long  as 
the  emission  does  not  swamp  out  the  directional  scattering 
effects. 


Alternate  Expression 

An  alternate  expression  for  obtaining  the  aerosol 
optical  thickness  from  sky  radiances  when  no  sun  scalar 
irradiance  measurement  is  available  can  be  developed. 
Substituting  Eq.  (4.25)  into  Eq.  (4.36)  and  rearranging  we 
get 

ytm(z)  “  4»r[  Li  (l,fl,,55  li) 

-  /.*  (j.8r  125/3)]  /  (  ,«<*)»!(•,)]  .  (4.37) 

The  sun  scalar  Irradiance  can  be  obtained  from  the 
average  sun  irradiance  out-of-the  atmosphere  ,«(°°)  by 
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,«<*)  -  J^(«)r»(x,o)'"(*l,  w?)J , 


where  ¥  is  the  angular  solar  diameter  for  that  date  and  'l' 
the  average  angular  solar  diameter.  Also  the  total 
tnnimltunce  can  be  expressed  as  a  function  of  the  com* 
ponent  optical  thicknesses  from  Eq,  (4.29) 

r»(3,0)-exp(  -/»/.(*)  -  w/„(x)  -  „/»(*)]  •  (4.39) 

Substituting  Eqs.  (4.38)  and  (4.39)  into  Eq.  (4.37) 
and  rearranging  we  get 

Wfw(*)expl  -MtK(i)m(8t)  ] 

-  4«r  Ui  (»,*,, JS  /J)  -  Lt.  (i.P,, 1250)1 
/[,f(«>)m(«1)('l'/ip)1exp{-m(4J)(  nim(i)  +  „/-(*)])]  .(4.40) 

Equation  (4.40)  cannot  be  solved  directly  for  the  aerosol 
optical  thickness  but  by  Interatlve  means  it  is  readily 
obtained.  It  assumes  a  reasonable  estimate  of  absorption 
optical  thickness  can  be  made  such  as  assuming  only 
oione  absorption  In  the  visible  or  that  a  value  can  be 
obtained  through  use  of  Program  LOWTRAN  at  other 
wavelengths. 

YtiKMon  SMIw 

The  two  methods  of  obtaining  the  total  optical  thick* 
nan  and  hence  transmittance  from  solar  almucantar  sky 
radiances  at  55  and  125  degrees  using  Eqa.  (4.37)  and 
(4.40)  with  Eq.  (4.29)  were  tested  using  the  1964  Visibility 
Laboratory  rooftop  data  previously  described  In  Section 
3.3. 

The  scanner  radiance  grid  included  almucantar 
measurements  within  ±0.7  degrees  for  all  four  Alters  for 
the  2  September  1964  nominal  70  degrees  sun  zenith 
angle  data  packages.  The  solar  transmiasometer  value  of 
transmittance  was  used  in  Eq.  (4.38)  to  obtain  the  value 
of  scalar  sun  irradiance  ,«(z)  for  use  In  Eq.  (4.37).  The 


values  of  and  basis  for  the  scalar  sun  irradiance  out  of  the 
atmosphere  ,<(°°).  the  Rayleigh  optical  thickness  and 
ozone  optical  thickness  are  given  in  Appendix  A. 

The  total  optical  thickness  /<«, (0)  and  total  transmit¬ 
tance  r-(0,0)  from  the  solar  transmissomcter  are  given 
for  each  Alter  in  Table  4.2.  The  ratio  of  the  optical  thick¬ 
ness  derived  using  Eq.  (4.37)  to  the  soiar  transmissometer 
value  is  given  in  column  4,  and  the  ratio  for  Eq.  (4.40)  in 
column  5.  Similarly,  the  ratio  of  the  derived  transmit- 
tances  to  measured  values  are  given  in  columns  7  and  8 
respectively.  The  values  compare  well  for  the  first  3  Alters 
and  less  well  for  Alter  4  meanwavelength  661  nanometers, 
although  even  there  the  transmittances  are  within  2  per¬ 
cent.  Livshits  and  Pavlov  (1970)  similarly  found  closer 
comparisons  at  the  shorter  wavelengths  for  mono¬ 
chromatic  measurements. 

Thus  the  shorter  method,  using  the  aimucantar  at 
55  and  125  degrees  with  or  without  an  Independent  sun 
Irradiance  measurement,  appears  to  be  valid  for  narrow 
band  and  broad  band  senaora  in  the  visible  spectrum. 

The  alternate  equation  Eq.  (4.40)  is  particularly  use¬ 
ful  with  airborne  scanner  data  where  independent  meas¬ 
urements  of  sun  Irradiance  are  not  available.  This  method 
was  used  to  obtain  the  aerosol  optical  thickness  for  flight 
C-466  airborne  data  at  6  and  0.2  kilometers  for  use  in  the 
example  in  Section  2.4  Tables  2.5  and  2.6. 

Since  this  is  an  approximative  method,  the  resultant 
aerosol  optical  thickness  precision  may  not  be  high,  but 
the  resultant  transmittances  have  high  precision  as  can  be 
seen  In  Table  4.2.  The  method  requires  an  unobscured 
sun  and  cloud-free  sky  at  55  and  125  degrees  from  the 
sun  in  the  almucantar. 


5.  SUMMARY 

Some  implications  of  the  equation  of  transfer  as  it 
relates  to  a  scattering  and  absorbing  medium  have  been 
explored.  The  m^jor  Implication  Is  that  a  measurement  of 
the  4ir  radiance  distribution  can  yield  a  great  deal  of  infor¬ 
mation  about  the  atmosphere  both  In  the  visible  and  the 
near  infrared.  If  the  solar  almucantar  Is  cloud-free  the 


Tabis  4.2,  Comparison  of  transmittance  determination  procedures 
(solar  slumcantsr  radiances  vs  xolsr  trensmlssometer). 
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aerosol  optical  thickness  can  be  derived.  1C  the  An  radi¬ 
ance  distribution  is  measured  at  several  altitudes,  the  radi¬ 
ance  arrays  can  be  tested  for  consistency  and,  if  con¬ 
sistent,  a  measure  of  absorption  obtained.  The  Visibility 
Laboratory  data  catalog  of  radiances,  both  airborne  and 
ground-based  are  a  mine  of  information  waiting  to  be 
tapped.  This  catalog  is  useful  for  testing  methods  of 
retrieval  of  basic  scattering  and  absorption  Information 
about  the  atmosphere  as  well  as  for  development  of 
atmospheric  models. 
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APPENDIX  A 


OPTICAL  AND  RADIOMETRIC 
CHARACTERISTICS  OF  SELECTED 
BROAD  BAND  SENSORS 

The  broad  band  sensors  referred  to  in  the  text  are 
muitipier  phototubes  with  spectral  Alters  interposed 
between  the  9ensor  and  the  optical  signal.  The  combina¬ 
tion  of  the  sensor  sensitivity  Sk  and  the  Alter  transmit¬ 
tance  Tk  is  the  resultant  sensitivity  of  the  Altered  photo¬ 
tube  SkTk.  The  standard  responses  which  each  optical 
system  attempts  to  duplicate  are  indicated  as  SkTk.  These 
relative  spectral  response  values  are  normalized  to  the 
peak  value. 


passband  8 X.  Similarly  the  radiometric  units  for  radiance 
of  W/sr  m1p.m  are  obtained  from  units  of  W/srm 7  by 
dividing  by  the  effective  passband. 

Inherent  Sun  Properties 

The  broad  band  values  of  sun  irradiance  out-of-the- 
atmosphere  at  mean  solar  distance  ,«(»)  are  obtained 
from  the  spectral  values  of  Johnson  (1954)  ,«(oo,x)  by 


1100 

1?(oo)  -  J?(oo,x)  Sktk  AX  /  8X  .  (A. 3) 


AVIZ  Filters 

The  relative  spectral  response  values  for  the  AVIZ 
Filters  2  through  6  are  graphed  in  Fig.  A-l  and  given  in 
tabular  form  in  Table  A.l.  Filter  code  9  is  the  true  photo- 
pic  response;  it  Is  included  in  Table  A.l  for  comparative 
purposes  only.  The  AVIZ  Aiters  are  applicable  to  airborne 
and  ground-based  data  measured  by  the  atmospheric  visi¬ 
bility  branch  of  the  Visibility  Laboratory  between  1970  and 
1978.  The  multiplier  phototubes  used  had  an  S-20  spec¬ 
tral  response. 


Fig.  A-l.  Standard  spectral  responses  for  AVIZ  (liters. 
Peak  Wavelengths  are:  2  -  473nm  Blue,  3  —  660nm  Red, 
4«S50nm  Photopic,  5-7S0nm  N.I.R.,  6«*440nm  S-20. 


The  spectral  average  inherent  sun  radiance  ,  r0  (X) 
was  computed  from 


sl0i\)  -,«(», x)  /  ns , 


(A.4) 


using  the  mean  angular  subtense  H,  of  6.819E-S  steradl- 
ans  based  on  a  mean  solar  distance  of  1.495E8  kilometers 
and  a  mean  solar  diameter  of  1.393E6  kilometers.  The 
spectral  center  sun  radiances  ,L„(X)  were  computed  using 
the  spectral  limb  darkening  factor  D(k)  based  on  the  limb 
darkening  functions  of  Mlnnaert  (1953) 


SL0M 


1 , Z0 (x)  /  DM  , 


(A, 5) 


Equations  similar  to  Eq.  (A. 3)  were  then  used  to 
obtain  the  broad  band  Inherent  sun  radiances  from  the 
spectral  values. 


Rayleigh  Atmosphere  Properties 

The  spectral  values  of  the  Rayleigh  total  volume 
scattering  coefficient  Ks  for  sea  level,  15  degrees  Celsius 
are  computed  from  (Eq.  14  from  Penndorf  (1957)), 

*j(X)  -  8ir3(nJ  -  1)J  (6  +  3 pn)  /  [  3x%(6  -  7pn ))  .(A.6) 


Spectral  Characteristics 

A  summary  of  the  response  characteristics  of  (he 
AVIZ  Alters  is  given  in  Table  A.2.  The  mean  wavelength 
X  is  deAned  by 

1100 

A-  £  X  AX/ 8  X  .  (A.l) 

The  effective  pAssband  8X  is  deAned  by 

1100 

8X  «■  L  $ktk  AX  .  (A. 2) 

inn  "  " 


The  number  density  t?  at  atmospheric  temperature  15 
degrees  Celsius  is  2.54743E19  cm"1  and  the  depolarization 
factor  p„  is  0.035.  The  refractive  index  n  can  be 
expressed  as  N’  where 

N'-n-l.  (A. 7) 

Now  rewriting  Eq.  (A.6)  in  terms  of  N'  and  with 
wavelength  in  nanometers 

*s(x)  -  3.443£12  w;4  nr1  (N'1  +  2 N')1  /  X4  .  (A. 8) 


The  radiometric  units  for  irradiance  of  are 

obtained  from  units  of  W/m 1  by  dividing  by  the  effective 
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The  refractive  term  /V'  can  be  expressed  as  a  function  of 
the  refractive  modulus  /V(0)  as 
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Tsble.A.l.  Relative  spectral  response  of  standards  for  AVIZ  Alters. 


FILTER  NUMBERS  AND  MEAN  WAVELENOTH 
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510 

* 

0.4140 

0.7930 

0.7100 

515 

♦ 

0.1)45 

0.7713 

0.7932 

550 

O.I4JO 

0.7500 

0.1610 

SJJ 

0.9045 

0.7130 

0.9149 

540 

0.9700 

0.7000 

0.9540 

545 

0.91)0 

0.4713 

0.9103 

550 

1.0000 

0.4570 

0.9950 

555 

0.944) 

0.4)15 

1.0001 

540 

0.91)0 

0.4100 

0.9930 

545 

4 

0.1415 

0.40)0 

0.9716 

570 

0.1040 

0.5140 

0.9530 

575 

0.719) 

0.5700 

0.9154 

510 

0.4)50 

0.5340 

01700 

515 

0.5515 

05115 

0.1161 

5f0 

0.4700 

0.3130 

0.7370 

595 

0.3950 

0.3040 

0.6949 

400 

0.1100 

0.4190 

0.6110 

405 

0.1410 

0.4750 

0.3641 

tto 

0.1040 

0.4410 

0.50)0 

FILTER  NUMBERS  AND  MEAN  WAVELENGTH 


No  2 

No  3 

No  4 

No  5 

No  6 

No  9 

Wive. 

Pseudo- 

True 

lenjth 

Blue 

Red 

Photcpic 

NIR 

S-20 

Phoiopic 

(nm) 

471  (nm) 

464  (nm) 

557  (nm) 

765  (nm) 

3)2  (nm) 

560  (nm) 

615 

. 

0.1610 

. 

0.4500 

0.4412 

610 

0.1300 

0.4390 

0.3110 

625 

4 

0.1035 

. 

0.4140 

0.3110 

630 

0.0000 

0.0110 

. 

0.41)0 

0.2630 

613 

0.0010 

0.0657 

4 

0.19)5 

0.2170 

640 

0.0416 

0.0504 

0  3740 

0.17)0 

645 

0.1791 

0.0411 

• 

0)343 

0.1)12 

650 

0.55)1 

0.0111 

0)150 

01070 

655 

0.9941 

0.0161 

0.3190 

0.0116 

640 

1.0000 

0.0111 

0)0)0 

0.0410 

445 

0.9421 

0.01  II 

. 

0.2143 

0.0446 

670 

0.1615 

0.0137 

0.2440 

0.0)20 

673 

0.7412 

0.01)9 

02410 

0.03)2 

610 

0.4774 

0.0120 

4 

0.2)00 

0.0170 

61) 

0.1515 

0.0105 

4 

0.1103 

0.0119 

690 

0.0495 

0.0090 

4 

0.1910 

0.0012 

69) 

. 

0.0146 

O.OOM 

4 

0.1735 

0.0057 

7(0 

0.0000 

0.0070 

0.1600 

0.0041 

70S 

* 

0.0041 

♦ 

0.1445 

0.0029 

710 

0.0051 

4 

0.1290 

0.0021 

71) 

0.0041 

0.1170 

0.0015 

710 

0.0041 

0.0000 

0.1050 

0.0010 

723 

0.00)1 

0.1005 

0.0911 

0.0001 

710 

0.001) 

0.2010 

0.0126 

0.000) 

735 

0.00)0 

0.4155 

0012) 

0.0004 

740 

0.0026 

06)00 

0.0619 

0.0001 

743 

0.0023 

0.1130 

0.0531 

0.0002 

730 

0.002) 

10000 

0.0497 

0.0001 

755 

0.0020 

0  9595 

0.0416 

0.0001 

740 

0.0011 

0.9190 

0.01)5 

0.0001 

765 

0.0017 

0.149) 

0.0292 

0.0000 

770 

0.0016 

0.7100 

0.0249 

77) 

0.0014 

0.6620 

0.0206 

710 

0.001) 

0.3440 

0.0162 

715 

0.0012 

0.4190 

0.0144 

790 

0.0012 

0.4)40 

0.0125 

795 

0.0011 

0.1720 

0.0107 

too 

0.0011 

0.3100 

0.0011 

105 

0  0005 

0.2673 

0.007) 

110 

0.0000 

0  2230 

0.0062 

113 

01123 

0.00)1 

120 

Table  A.2.  Spectral  characteristics  summary  for  AVIZ  Alters. 


Spectral  Cturactarlillci 

Inhirani  Sun  Propaniei  IJohnion  (19)41) 

Rtyiel|h  Aimoephtra  Fropenlei  (13*0 

Otont 

Optical 

Thlcknaii 

a 

Pillar 

Cork 

No. 

P*ik 

Wivelanuh 

(nm) 

Maen 

Wivaitniih 

(nm) 

EITactlva 

Paaeband 

(nm) 

trradlanca 

(M'/»,m«i) 

Rid  line*  ( 

Auanuatlon 

Lansih 

(m) 

Total  Sctiieilnt 
CoalActam 

(per  mi 

Vertical 

Reditnce 

Trtnemiiitnct 

Averts* 

Canter 

2 

475 

411 

19.9 

2.I4E+03 

3  1  IE +07 

4.01E+07 

4.I4E+04 

2  07E-05 

0139 

3.71E-J 

3 

460 

664 

102 

I37E+03 

2.30E+O7 

2.7SE+07 

I.B6E+05 

3  4IE-06 

0955 

1  SILO 

4 

350 

357 

115 

1  90E+03 

27IE+07 

147E+07 

I.91E+04 

1  I5E-05 

0  907 

)  OIE-2 

3 

750 

763 

50.4 

123  E  -r  03 

1 10E+07 

2.I0E+07 

3.3IE+OS 

1 OIE-06 

0  974 

j  1  61E-3 

6 

440 

3)2 

111.3 

1  9IE+03 

2I0E+07 

1.33E+07 

7.22E+04 

1  64E-03 

n  167 

1 

9 

353 

340 

106.9 

I.I9E+03 

2.77E+07 

3.4JE+07 

9.22E+04 

1  I3E-03 

0907 

1  3  07E-I 

L _ 

Ls'.oAv 


•32- 


i_v *j\.kT 


APPENDIX  A 


/V'-/Y(0)l£-6.  (A. 9) 

The  refractive  modulus  was  computed  using  the  dispersion 
formula  of  Edlen  (1953)  for  the  refractive  modulus  N(0) 
at  15  degrees  Celsius  at  sea  level 

N(0)  -  64.328  +  29498.10  /  (146  -  1  /  **) 

+  255.40/  (41  -1/AJ),  (A.  10) 


where  X  is  in  pm.  The  spectral  Rayleigh  vertical  beam 
transmittance  through  the  atmosphere  may  be  computed 
by  means  of  the  scale  height  H 


aT^M-expl-jito)//) .  (A.11) 


The  scale  height  for  the  V.  S.  Standard  Atmosphere  (1976) 
which  has  a  sea  level  temperature  of  15  degrees  Celsius,  is 
8434.7  meters. 

Attenuation  Length.  The  attenuation  length  £  is 
defined  as  that  distance  at  which  the  signal  is  attenuated  to 
1/e.  So  monochromatically,  the  Rayleigh  attenuation 
length  is 

*£(A)-l/fti<X>.  (A.12) 

A  broad  band  attenuation  length  is  similarly  defined  as  the 
distance  at  which  the  signal  is  attenuated  to  1/e.  Thus 

exp(-l)-exp(~r/j(£) 

1100 

■  ^  axpl-^sfA)  J  SJ^Ak/hk  .  (A.  13) 


Equation  (A.  13)  cannot  be  eolved  directly  for  the  Rayleigh 
attenuation  lcr.jjth  *£  since  the  range  r  is  both  unknown 
and  in  the  defining  term.  However,  it  can  be  solved  by 
Iteration  within  acceptable  limits. 

Volume  Scattering  Coefficient.  The  Rayleigh  volume 
scattering  coefficient  was  needed  as  a  lower  limit  for  com¬ 
parison  to  the  values  measured  by  the  integrating 
nephelometer.  Therefore,  the  limiting  value  a  5  is  defined 

as 


1100 

^  Ms(k)MiK55WtCSJZ*k 
1100 

/  J  AfU.SSOO  K)  5^77  Ak  .  (A, 14) 


A  Xenon  source  is  used  in  the  nephelometer,  therefore 
the  radiant  exitance  (omittance)  in  Eq.  (A.  14)  was  approx¬ 
imated  by  a  blackbody  radiator  of  5500  degrees  Kelvin 
Af(k,3500'K).  The  radiant  exitance  M(k,T)  is 


At{k,T)  m  itHk,T)  .  (A.15) 


The  Hk.T)  is  computed  from  the  classical  equation 
(Wolfe  (1978)  Table  1.7) 


t(A.r)-r,  /  |irA5(c-*-  1)1  .  (A. 16) 

where 


r,  -  2« Ac3  -  3.741382£-16  Wm 3  (A,  1 7) 


and 

x-cj/UT).  (A.18) 


The  c  is  the  speed  of  light,  h  is  the  Planck  constant  and 
Cj  -  Ch  /  k  -  1.438786£-2  mK  (A.19) 


where  k  is  the  Boltzmann  constant.  The  constants  are 
from  Driscoll  and  Vaughn  (1978)  Table  A.l. 

Transmittance.  The  Rayleigh  transmittance  through 
the  atmosphere  was  needed  as  a  limit  for  comparison  to 
the  transmittance  measurements  made  with  a  solar 
transmissometer.  Therefore  the  broad  band  Rayleigh 
transmittance  was  defined  as, 

uoo 

n  r.(0.0)  *  £  ,L0(k)  h  TJk)  ffX  Ak 

UOO 

/  J  ,Ltlk)S^  Ak  .  (A, 20) 

where  ,Mx)  is  the  spectral  inherent  center  sun  radiance 
as  defined  by  Eq.  (A.5). 


Ozone  Absorption 

The  spectral  values  of  ozone  absorption  a  per  cm 
for  the  Chappuis  bands  are  taken  from  the  tabular  values 
of  Vigroux  (1953)  which  are  in  good  agreement  with  the 
Inn  and  Tanaka  (1953)  values.  Since  the  ozone  absorp¬ 
tion  was  to  be  used  in  conjunction  with  the  solar  transmis¬ 
someter  values,  the  absorption  a  per  cm  for  each  filter 
was  computed  from  the  spectral  values  by 

uoo 

0“  5)  |,(“A)  4 A 

MOO 

/  ^  ,i(°®,A)  Tk  Ak  .  (A. 21) 

The  total  ozone  was  from  the  U.  S.  Standard  Atmo¬ 
sphere  (1976)  which  was  0.345  atm-cm.  The  ozone  opti¬ 
cal  thickness  was  computed  for  each  filter  from 

r^(0.0>  -  exp  t  -  ar„(0)l  .  (A.22) 
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Rooftop  Filters 

The  relative  spectral  response  values  for  the  broad 
band  sensors  used  on  the  rooftop  of  the  Visibility  Labora¬ 
tory  buildini  during  1964  are  graphed  in  Fig.  A-2  and 
given  in  tabular  form  in  Table  A.3.  The  multiplier  photo¬ 


tubes  had  an  S-ll  spectral  response.  A  summary  of  the 
response  characteristics  of  the  Rooftop  1964  filters  is  given 
in  Table  A.4. 

Rooftop  Filter  3  is  the  photopic  sensor.  The  photo¬ 
tubes  were  carefully  filtered  so  that  they  closely  approxi¬ 
mated  the  photopic  sensor. 


Table  A.3.  Relative  spectral  response  of  standards  for  rooftop  1964  filters. 


FILTER  NUMBERS  ANO  MEAN  WAVELBNOTH 


FILTER  NUMBERS  AND  MBAN  WAVELBNOTH 


Wiv«t*n(th 

(nm) 


W»y»kfl|lh 

(nm) 

No.  1 

459  (nm) 

No.} 

505  (nm) 

No.  3 

560  (nm) 

NO.  4 

661  (nm) 
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000 w 

04405 

0.1161 

SK> 

0.0027 

0.4130 

0.7544 

. 

$95 

0.0000 

0.3420 

0.6941 

. 

600 
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0.6309 

0.0000 

60S 

03940 

0.5647 

0.0110 

610 

0.1970 

0.5039 

0.1740 

61 S 

01390 

0.4411 

0.5030 

630 

0.1310 

0.3109 

0.1310 

US 

o  toot 

0.3309 

0.9140 

6  JO 

0.0791 

0.3649 

1.0000 

US 

0.06R4 

0.3170 

0.9410 
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00376 

0.1750 
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00491 
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0.797J 

6S0 

0  0430 
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0.0JJ4 

0.0116 
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0.0610 
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66S 
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0.0191 
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US 
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001 1 1 
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699 

oom 

0.0057 
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001  to 

0.0041 

0.1110 

70S 

0.0141 

0.0039 

0.3500 

710 

00173 

00031 

0.3730 
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. 

0.01  SI 

0.0013 

0.3355 

730 

0014] 

0  0010 

01990 

73S 

00144 

0.0007 

0.1195 

730 

0.0144 

0.0005 

o.iaoo 

7J5 

0.0141 

0.0004 

01650 

740 

0.0135 

00003 

0.1440 
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00133 

0.0003 

0.1363 

7S0 

0.0109 

0.0001 

0.1070 

7S5 

0.0094 

0.0001 

0.0913 

760 

0.0013 
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0.0754 

765 

0.0061 

ooooo 

0  0421 

770 

0.0055 

0.0301 

775 
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00405 

710 

00036 
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715 
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00359 
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00030 

0.0117 

795 

0  001 5 

00137 
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0.001 1 

0  0016 
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n«.  A-2.  SUndtrd  spectral  reaponaea  for  Visibility  Laborttory 
rooftop  Alters. 


TabU  A.4.  Spectral  characteristica  summary  for  rooftop  1964  niters. 
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GLOSSARY  AND  NOTATION 


The  notation  used  in  reports  and  journal  articles  produced  by  the  Visibility  Laboratory  staff  follows,  in  general,  the 
rules  set  forth  in  pages  499  and  500,  Duntiey  ti  al,  (1957).  These  rules  are: 

•  Each  optical  property  is  indicated  by  a  basic  (parent)  symbol. 

•  A  presubecript  may  be  used  with  the  parent  symbol  as  an  identifier,  e,g,,  b  indicates  background  while  t  denotes 
an  object,  /.#.  target, 

•  A  poetsubecript  may  be  used  to  indicate  the  length  of  a  path  of  sight,  r  denotes  an  apparent  property  as 
meaiutod  at  the  end  of  a  path  of  sight  of  length  r,  while  o  denotes  an  inherent  property  based  on  the  hypothetical 
concept  of  a  photometer  located  at  zero  distance  from  an  object,  /,#.  target. 

•  A  postsuperacript*  or  postsubscript*,  is  employed  as  a  mnemonic  symbol  signifying  that  the  radiometric  quantity 
has  been  generated  within  the  path  or  path  segment  by  the  scattering  of  ambient  light  reaching  the  path  from  all 
directions  and/or  by  emission. 

•  The  parenthetical  attachments  to  the  parent  symbol  denote  altitude  and  direction.  The  letter  z  indicates  altitude  in 
general;  z,  is  used  to  specify  the  altitude  of  a  target.  The  direction  of  a  path  of  sight  is  specified  by  the  zenith 
angle  9  and  the  azimuth  In  the  case  of  irradiances,  the  downwelling  irradlance  is  designated  by  </,  the 
upwelllpg  by  u, 

•  The  radiometric  symbols  used  herein  now  correspond  to  the  OSA  recommendations  in  Section  1  of  Driscoll  and 
Vaughn  (1971).  Prior  to  June  1980,  the  symbol  used  for  radiance  L  was  /V,  for  irradlance  £  was  //,  and  for 
attenuation  length  fwas  L, 


Symbol 


Quantity 


Albedo  at  altitude  z,  defined 
«4(i)  -  £(z,w)  /  E(t,d). 


Absorption  coefficient. 


sr~]  m“' 


Diffuse  component  of  the  solar  almucantar  sky  radiance  ratio  to  the  sun  scalar 
irradlance  and  relative  alrmass 

C(z)  -  f'Ux,9\4')(r{ztp')d(\  /  ,*(z)  +  L,a(t)  /  ,«(z). 


Speed  of  light  c  -  2.99792458£8. 


First  radiation  constant  o  -  Inhv*  -  3.74I832E-16. 


Second  radiation  constant  cj  -  ch/'k  -  1.438786£-2. 


Radiance  distribution  function  D(z)  -  ^L(t,9\(b')sece'dCl  /  «(z) 

Limb  darkening  factor  relating  tho  average  sun  radiance  and  the  center  sun 
radiance  D(k)  -  ,L0  /  ,L0. 

Spectral  irradlance  (formerly  symbol  H)  defined  as  £„  -  j L^2,9,<b)  cos 9  di\. 
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Symbol  Units 


Quantity 


E 

Broad  band  sensor  irradiance,  defined  as  £  -  j  £»  5*  7\  dk  /  6k. 

0 

E(x,d) 

W/m 1 

Irradiance  produced  by  downwelling  flux  as  determined  on  a  horizontal  flat  plate 
at  altitude  z  [formerly  H(x,d)].  In  this  report  d  is  used  In  place  of  the  minus 
sign  in  the  notation  !//(*,—)]  which  appears  in  Duntley  (1969).  This  property 
may  be  defined  by  the  equation  E(x,d)  -  j  L(x, 9', d>')  cosP'dn’. 

E(x,u) 

W/m 1 

Irradiance  produced  by  the  upwelllng  flux  as  determined  on  a  horizontal  flat  plane 
at  altitude  z  (formerly  //(z,m)].  Here  u  Is  substituted  for  the  plus  sign  formerly 
used  in  the  notation  lff(x,+)l 

Hit) 

m 

Scale  height  at  altitude  z,  the  height  of  a  homogeneous  atmosphere  having  the 
density  of  the  layer  at  altitude  x. 

A 

J/s 

Planck  constant  A  -  6,6261 76£-34, 

k 

J/K 

Boltzmann  constant  k  -  1.380662E-23. 

Li 

W/v  m* 

Spectral  radiance  (former  symbol  N). 

L 

BVzr 

Broad  band  sensor  radiance  is  defined  as  L  -  J  L*  Skfk  dk  /  6k. 

• 

Lj(  z<,0,d) 

H7sr  m1 

Inherent  radiance  based  on  the  hypothetical  concept  of  a  photometer  located  at 
zero  distance  from  an  object  at  altitude  z  in  the  direction  specified  by  zenith  angle 
9  and  azimuth  d- 

W/v  m1 

Apparent  radiance  as  determined  at  altitude  z,  from  the  end  of  a  path  of  sight  of 
length  r  at  zenith  angle  9  and  azimuth  *.  This  property  may  be  defined  by 
Lr(s,9,d>)  -  t,(z„M)r,(z,g)  +  l;u, M). 

L(\,T) 

M7»  «* 

Black  body  radiance  at  wavelength  k  and  temperature  T. 

IP/ir  w* 

Equilibrium  radiance  at  altitude  z  with  the  direction  of  the  path  of  sight  specified 
by  zenith  angle  9  and  azimuth  d<  This  property  is  a  point  function  of  position 
and  direction. 

As  discussed  by  Duntley  *i  al.  (1957),  many  image  transmission  phenomena  are 
most  clearly  understood  in  terms  of  the  concept  of  equilibrium  radiance.  This 
concept  is  a  natural  consequence  of  the  equation  of  transfer  which  states 
analytically  that  in  any  path  segment  the  difference  between  the  output  and  Input 
radiances  is  attributable  to  a  gain  term  and  a  loss  term,  such  that  some  unique 
equilibrium  radiance  £f(z,M)  must  exist  at  each  point  such  that  the  loss  of 
radiance  within  the  path  segment  Is  exactly  balanced  by  the  gain, 
A/.(z,0,d)  -  0. 
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Symbol  Units 


Quantity 


L,(tt9,4)  W/v  ml 

Lsto,9A)  W/sr  m> 


By  virtue  of  this  concept  and  the  equation  of  transfer,  one  can  show  that  each 
segment  of  every  path  of  sight  has  associated  with  it  an  equilibrium  radiance,  and 
that  the  space  rate  of  change  in  image  forming  radiance  caused  by  the  path 
segment  is  in  such  a  direction  as  to  cause  the  output  radiance  to  be  closer  to  the 
equilibrium  radiance  than  is  the  input  radiance.  This  segment  by  segment 
convergence  of  the  apparent  radiance  of  the  object  field  to  the  dynamic 
equilibrium  radiance  was  clearly  illustrated  by  the  data  in  the  1957  paper 
referenced  above. 


Path  function  at  altitude  z  with  the  direction  of  the  path  of  sight  specified  by 
zenith  angle  9  and  azimuth  4.  This  property  is  defined  by  the  equation 
+  Ua(z). 

Path  function  due  to  scattering,  defined  by  L,(z,M)  -  f  orh,p')Lh,»\4')d(\'. 


L*<i) 


W/sr  mi  Path  function  due  to emittance  (absorption)  L,fl(z)  -  a(z)  L(K,T). 


L/(r,0,#)  W/sr  mJ 


Path  radiance  as  determined  at  altitude  z  at  the  end  of  a  path  of  sight  of  length  r 
in  the  direction  of  zenith  angle  9  and  azimuth  4- 


LL(z,9,4)  W/sr  m 1 


Sky  radiance  at  altitude  1,  zenith  angle  9  and  azimuth  4-  Also  the  path  radiance 
for  the  path  of  sight  of  length  «>  from  out  of  the  atmosphere  to  altitude  z. 


,Lm(zt9 ,,4)  W/sr  ml  Apparent  radiance  of  the  center  of  the  solar  disk  u  determined  from  the  end  or  a 

path  of  sight  of  length  «°  from  out  of  the  atmosphere  to  altitude  z  at  the  zenith 
angle  of  the  sun  9,. 


S(s)  m 


Attenuation  length  at  altitude  z.  The  attenuation  length  is  the  distance  at  which 
the  signal  is  attenuated  to  lfe. 


W/m 1 

m(Kt) 

W/m * 

M«(s) 

W/m ‘ 

»«>(*) 

W/m 1 

m^(z) 

W/m> 

mn(i,«) 

kg/m 1 

Black  body  exitance  (emittance)  AfU.D  *  irZ,(\,D. 

Black  body  scalar  exitance  (emittance)  m(K,  T)  -  4ir L(k,  T). 

Scalar  exitance  (emittance)  per  length  m,(z)  «  £L,(z,9,4)d 0, 

Scalar  exitance  (emittance)  per  length  due  to  scattering 
*U((z)  «  ^L^{z,9,4)d(\  -«(z)s(z), 

Scalar  exitance  (emittance)  per  length  due  to  absorption  m,g  (z)  -  a  (r)m(X,  T). 
Absolute  air  mass  at  angle  9  mm(z,9)  -  J  p(t)dr. 
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Symbol 

Units 

Quantity 

mM(z,P)/m..(z,0^ 

none 

Relative  optical  air  mass. 

mii,) 

nom 

Relative  optical  airmass  at  zenith  anile  of  the  sun,  shorthand  for 
/*»•(*,  9,)/mmit%  0*9, 

Nit) 

none 

Refractive  modulus  Nit)  -  [n(z)  -  l)£6. 

Hit ) 

none 

Refractive  index. 

Oil) 

none 

Optical  scattering  mixing  ratio  at  altitude  t.  This  quantity  la  defined  as  the  ratio 
of  the  total  volume  scattering  coefficient  at  altitude  t,  to  the  molecular  (or 
Rayleigh)  volume  scattering  coefficient  at  the  same  altitude  t. 
Qit)  -  tit)/ *  sit). 

aM) 

none 

Volume  scattering  function  ratio  at  altitude  t.  This  quantity  is  defined  as  the  ratio 
of  the  total  volume  scattering  function  at  altitude  s  and  scattering  angle  0,  to  the 
molecular  (or  Rayleigh)  volume  scattering  function  at  the  same  altitude  and 
scattering  angle.  9(z,0)  -  <r(z,0)  /  *<r(z,0). 

r 

m 

Path  length,  for  paths  of  sight  at  zenith  angles  0  to  70  degrees,  r  -  step  Az. 

sx 

none 

Standardized  relative  spectral  response  of  fllter/cathode  combination  where  Is 
spectral  sensitivity  of  the  multiplier  phototube  cathode  and  7\  is  spectral 
transmittance  of  optical  Alter.  The  relative  spectral  response  values  are 
normalized  to  the  peak  value. 

i(m) 

m~l 

Total  volume  scattering  coefficient  as  determined  at  altitude  s.  This  property  may 
be  deAned  by  the  equations 

i(z)  -  £o(x,0)rfO  -  *s(z)  +  ntsit). 

In  the  absenoe  of  atmospheric  absorption,  the  total  volume  scattering  coefficient  is 
numerically  equal  to  the  attenuation  coefficient. 

mil) 

m~' 

Volume  scattering  coefficient  for  Mie  it.  particulate  or  droplet,  scattering  at 
altitude  z. 

mil) 

m-' 

Volume  scattering  coefficient  for  Rayleigh  /.#.  molecular,  scattering  at  altitude  z. 

T 

K 

Absolute  temperature  in  degrees  Kelvin. 

none 

Radiance  transmittance  as  determined  at  altitude  z  for  a  path  of  sight  of  length  r 
at  zenith  angle  0  (formerly  referred  to  as  "beam”  transmittance).  This  property  is 
independent  of  azimuth  in  atmospheres  having  horizontal  uniformity.  It  is  always 
the  same  for  the  designated  path  of  sight  or  Its  reciprocal. 

none 

Radiance  transmitancc  for  the  path  of  sight  at  zenith  angle  B  from  out  of  the 

atmosphere  to  the  altitude  z. 


•39- 


APPENDIX  B 

Symbol 

Units 

mi 

l| 

Quantity  ' 

'**<*) 

none 

Optical  thickness  /41(z)  -  J*  a{z)dz. 

'»  ; 

tJ 

•  *S 

none 

ii 

Absorption  optical  thickness  fl/A,(z)  -  f  a(z)<tz.  • 

•i  ' 

*/*,(«) 

none 

] 

Mie  optical  thickness  «rAI(z)  «  /  «z(z)<fc.  * 

*t  t 

; 

*'*»<*) 

none 

Rayleigh  or  molecular  optical  thickness  *tkl(z)  -  f  ns(z)dz , 

w 

J/k' 

1 

Radiant  density. 

* 

w 

W 

• 

t 

Watt  w  -  I/s.  ; 

S 

i 

m 

! 

Altitude,  usually  used  as  above  ground  level. 

• 

* 

m 

i. 

Altitude  of  any  applicable  object  or  target. 

«(») 

m"‘ 

t 

Volume  attenuation  coefficient  u  determined  at  altitude  z.  a(z)  -  a(z)  +  s(z).  ; 

t  * 

0 

d«| 

Symbol  for  scattering  angle  of  flux  from  a  light  source.  It  is  equal  to  the  angle 
between  the  line  from  the  source  to  any  unit  scattering  volume  and  the  path  of  a 
ray  scattered  off  this  direct  line.  See  illustration. 

,  l  S 

'yk/  # 

—— ^Sun j — 

j^Unit  Volume 

X  / 

/V  /  O 

yS  7  $ 

A 

none 

Symbol  to  Indicate  incremental  quantity  and  used  with  r  and  z  to  Indicate  small, 
discrete  increments  In  path  length  r  and  altitude  r. 

ax 

nm 

* 

Effective  passband  (formerly  designated  "response  area”)  for  a  Altered  sensor  is 

defined  as  8A  -  J  ST77  dk. 

0 

\V 
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Symbol 

Units 

Quantity 

€<z) 

W/m1 

Scalar  irradiance.  This  may  be  defined  as  the  radiant  flux  arriving  at  a  point,  from 
all  directions  about  that  point,  at  altitude  z  (Tyler  and  Preisendorfer  (1962)). 
«(z)  -  j'  L(z,9’,<y)dn', 

*«(i) 

W/m 1 

Diffuse  scalar  irradiance  rf«(z)  -  t(z)  -  ,«(z). 

!«(«) 

W/m 1 

Sun  scalar  irradiance  at  altitude  z  ,«(z)  -  ,«(<»)  r»(z,9,). 

( 

m 

Radius  of  the  earth. 

0 

deg 

Symbol  for  zenith  angle.  This  symbol  is  usually  used  as  one  of  two  coordinates  to 
specify  the  direction  of  a  path  of  sight. 

0' 

deg 

Symbol  for  zenith  angle  usually  used  as  one  of  two  coordinates  to  specify  the 
direction  of  a  discrete  portion  of  the  sky. 

•* 

deg 

Zenith  angle  of  the  sun. 

X 

nm 

Symbol  for  wavelength. 

j: 

nm 

Mean  wavelength  Is  defined  as  k  -  J  X  SkTk  dk  /  OX. 

• 

sr> 

Ratio  of  solar  almucantar  sky  radiance  to  scalar  sun  irradiance  and  airmass. 
m(a i/S)  -  9)/  ( jfGrimf#,)]. 

Kl) 

H7m* 

Net  Irradicnce.  ((z)  -  E(t,d)  -  £(z,u). 

p(z) 

kg/m* 

Density  at  altitude  z. 

/w'lr* 

Symbol  for  volume  scattering  function.  Parenthetical  symbols  are  z  to  designate 
altitude  and  /3  to  designate  the  scattering  angle  from  a  source. 

cr(i,/9)/j  (x) 

jr1 

Normalized  volume  scattering  function.  This  may  be  definod  by  the  equation 
^f(  <z(z,0)  /  j(z)]dfi  -  1. 

jr 

Optical  thickness  function.  ta,(z,0)-  f<r(z,0)dz 
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Symbol  Units  Quantity 


M  dog  Symbol  for  azimuth.  The  azimuth  is  the  angle  in  the  horizontal  plane  of  the 

observer  between  a  fixed  point  and  the  path  of  sight.  The  fixed  point  may  be  for 
example,  true  North,  the  bearing  of  the  sun.  or  the  bearing  of  the  moon.  This 
symbol  is  usually  used  as  one  of  two  coordinates  to  specify  the  direction  of  a  path 
of  sight. 


p*  deg  This  symbol  for  azimuth  is  usually  used  as  one  of  two  coordinates  to  specify  the 

direction  of  a  discrete  portion  of  the  sky. 


Mr  deg  Angular  solar  radius  at  true  earth-to-sun  distance. 


$  deg  Angular  tolar  radius  at  mean  solar  distance. 


0  sr  Symbol  for  solid  angle.  For  a  hemisphere,  H  -  2w  steradians.  For  a  sphere, 

n  -4tr  steradians. 


«(») 


none 


Single  scattering  albedo  u(r)  -  j(z)  /  a(z). 
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Table  B.l.  Notational  equivalencies. 


Chandrasekhar  (1960) 

Visibility  Laboratory 

Symbol 

Definition 

Symbol 

Definition 

E  .. 

Radiant  energy 

E 

Irradiance 

irf 

Net  flux 

{ 

Net  irradiance 

l 

Specific  intensity 

JLGr.M) 

Radiance 

J 

Average  intensity 

«/4ir 

Average  radiance 

Mass  emission  coefficient 

L.(x.M) 

J 

p(x) 

Path  function/density 

Jp 

Mass  emission  coefficient  x  density 

l.(x.M) 

Path  function 

jU) 

Mass  emission  coefficient  due 

Path  function  due  to  scattering 

to  scattering 

p(z) 

divided  by  the  density 

k 

Mass  absorption  coefficient 

alp 

Attenuation  coefficient  divided 
by  the  density 

kp 

Maas  absorption  coefficient  x  density 

a(i) 

Attenuation  coefficient 

p(ccs9) 

Phase  function 

4ir<r(fl)/a 

Four  n  times  volume  scattering 
function/attenuation  coefficient 

» 

Thickness 

r 

Path  length 

z 

Source  function 

Equilibrium  radiance 

u 

Integrated  energy  density 

w 

Radiant  density 

“o 

Single  scattering  albedo 

«(x) 

Single  scattering  albedo 

0 

Scattering  angle 

p 

Scattering  angle 

& 

Polar  angle 

9 

Zenith  angle 

M 

coeine  of  polar  angle 

cos  9 

cosine  of  zenith  angle 

a 

Frequency 

1/A 

Inverse  wavelength 

P 

Density 

P 

Density 

T 

Normal  optical  thickness 

'at 

Optical  thickness 

9 

Atlmuth 

9 

Azimuth 

« 

Solid  angle 

n 

Solid  angle 
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